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Methods for the design of physical parameterization schemes that possess certain invari- 
ance properties are discussed. These methods are based on different techniques of group 
classification and provide means to determine expressions for unclosed terms arising in the 
course of averaging of nonlinear differential equations. The demand that the averaged equa- 
ls! tion is invariant with respect to a subalgebra of the maximal Lie invariance algebra of the 
^ unaveraged equation leads to a problem of inverse group classification which is solved by 
Q the description of differential invariants of the selected subalgebra. Given no prescribed 
symmetry group, the direct group classification problem is relevant. Within this framework, 
the algebraic method or direct integration of determining equations for Lie symmetries can 
be applied. For cumbersome parameterizations, a preliminary group classification can be 
carried out. The methods presented are exemplified by parameterizing the eddy vorticity 
flux in the averaged vorticity equation. In particular, differential invariants of (infinite di- 
mensional) subalgebras of the maximal Lie invariance algebra of the unaveraged vorticity 
equation are computed. A hierarchy of normalized subclasses of generalized vorticity equa- 
tions is constructed. Invariant parameterizations possessing minimal symmetry extensions 
are described. A restricted class of invariant parameterization is exhaustively classified. The 
^ physical importance of the parameterizations designed is discussed. 

<N 

> 1 Introduction 

O 

The problem of parameterization is one of the most important issues in modern dynamic mete- 
C<1 orology and climate research [16, 42]. As even the most accurate present days numerical models 

are not capable to resolve all small scale features of the atmosphere, there is a necessity for find- 
^""i ing ways to incorporate these unresolved processes in terms of the resolved ones. This technique 

is referred to as parameterization. The physical processes being parameterized in numerical 
^ weather and climate prediction models can be quite different, including e.g. cumulus convection, 

• ^ momentum, heat and moisture fluxes, gravity wave drag and vegetation effects. The general 

^ problem of parameterization is intimately linked to the design of closure schemes for averaged (or 

5^ filtered) nonlinear equations. By averaging, a nonlinear differential equation becomes unclosed, 

that is, there arise additional terms for which no prognostic or diagnostic equation exist. These 

terms must hence be re-expressed in a physically reasonable way to be included in the averaged 

equations. 

It has been noted in [43] that every parameterization scheme ought to retain some basic 
properties of the unresolved terms, which must be expressed by the resolved quantities. These 
properties include, just to mention a few, correct dimensionality, tensorial properties, invariance 
under changes of the coordinate system and invariance with respect to Galilean transformations. 
While the formulation of a parameterization scheme with correct dimensions is in general a 
straightforward task, not all parameterization schemes that have been used in practice are 
indeed Galilean invariant. An example for this finding is given by the classical Kuo convection 
scheme [19, 20]. In this scheme, it is assumed that the vertically integrated time-change of 
the water vapor at a point locally balances a fraction of the observed precipitation rate [8, 
pp. 528]. This also implies that the moisture convergence is proportional to the precipitation 
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rate. However, while the precipitation rate is clearly a Galilean invariant quantity, the moisture 
convergence depends on the motion of the observer [17]. That is, the Kuo scheme does not 
properly account for pure symmetry constraints, which is a potential source of unphysical effects 
in the results of a numerical model integration. 

The latter finding is the main motivation for the present investigations. Galilean invariance 
is an important example for a Lie symmetry, but it is by no means the only invariance character- 
istic that might be of importance in the course of the parameterization process. This is why it 
is reasonable to focus on parameterization schemes that also preserve other symmetries. This is 
not an academic task. Almost all real-world processes exhibit miscellaneous symmetry charac- 
teristics. These characteristics are reflected in the symmetry properties of differential equations 
and correspondingly should also be reflected in case where these processes cannot be explicitly 
modeled by differential equations, i.e. in the course of parameterizations. What is hence desir- 
able is a constructive method for the design of symmetry- preserving parameterization schemes. 
It is the aim of this paper to demonstrate that techniques from group analysis do provide such 
constructive methods. In particular, we state the following proposition: 

Any problem of finding invariant parameterizations is a group classification problem. 

Implications following from the above proposition form the core of the present study. It appears 
that this issue was first opened in [29], dealing with the problem of turbulence closure of the 
averaged Navier-Stokes equations. We aim to build on this approach and extend it in several 
directions. As the equations of hydrodynamics and geophysical fluid dynamics usually possess 
wide symmetry groups [2, 4, 5, 10, 12], the design of symmetry-preserving parameterizations 
will in general lead to a great variety of different classes of invariant schemes. 

Needless to say that the parameterization problem is too comprehensive both in theory and 
applications to be treated exhaustively in a single paper. Therefore, it is crucial to restrict to a 
setting that allows to demonstrate the basic ideas of invariant parameterizations without overly 
complicating the presentation by physical or technical details. This is the reason for illustrating 
the invariant parameterization procedure with the rather elementary barotropic vorticity equa- 
tion. For the sake of simplicity, we moreover solely focus on local closure schemes in the present 
study. That is, the quantities to be parameterized at each point are substituted with known 
quantities defined at the same respective point [43]. This renders it possible to thoroughly use 
differential equations and hence it will not be necessary to pass to integro-differential equations, 
as would be the case for nonlocal closure schemes. On the other hand, this restriction at once 
excludes a number of processes with essential nonlocal nature, such as e.g. atmospheric con- 
vection. Nevertheless, there are several processes that can be adequately described within the 
framework of the present paper, most notably different kinds of turbulent transport phenom- 
ena. 

The organization of this paper is the following: Section 2 discusses different possibilities for 
the usage of symmetry methods in the parameterization procedure, most noteworthy the appli- 
cation of techniques of direct and inverse group classifications. We restate some basic results 
from the theory of group classification and relate them to the parameterization problem. Sec- 
tion 3 is devoted to the construction of several parameterization schemes for the eddy vorticity 
flux of the vorticity equation using the methods introduced in the previous section. Bases of 
differential invariants and operators of invariant differentiation for subalgebras of the maximal 
Lie invariance algebra of the vorticity equation are computed and used in the framework of 
invariant parameterization (Section 3.1). It should emphasized that up to now only very few 
examples on exhaustive descriptions of differential invariants for infinite dimensional Lie alge- 
bras exist in the literature [6, 11]. A hierarchy of nested normalized subclasses of a class of 
generalized vorticity equations is constructed in Section 3.3. Additionally, in Section 3.2 the 
equivalence algebras of some subclasses are directly found within the framework of the infinites- 
imal approach. The algebraic method of group classification is used to determine inequivalent 
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invariant parameterization schemes. For a restricted class of generalized vorticity equations, 
it is proved in Section 3.4 that the algebraic method provides an exhaustive description of all 
inequivalent parameterizations of the eddy vorticity flux. For a wider class of generalized vor- 
ticity equations, in Section 3.5 we study the problem of invariant parameterization within the 
framework of preliminary group classification. Namely, inequivalent invariant parameterizations 
possessing at least one-dimensional symmetry extensions are listed. A short discussion of the 
results of the paper is presented in Section 4, together with an outlook on forthcoming works 
in the field of invariant parameterization theory. In Appendix A, details on the classification of 
inequivalent one-dimensional subalgebras of the equivalence algebra from Theorem 1, which is 
used in Section 3.5, can be found. 

2 The general idea 

Throughout the paper, the notation we adopt follows closely that presented in the textbook [31]. 
Let there be given a system of differential equations 

A'(rE,n(„)) = 0, l = l,...,m, (1) 

where x = {x^, . . . , x^) denote the independent variables and the tuple includes all dependent 
variables u = (u^, . . . , u'^) as well as all derivatives of u with respect to x up to order n. Hereafter, 
subscripts of functions denote differentiation with respect to the corresponding variables. 

Both numerical representations of (1) as well as real-time measurements are not able to 
capture the instantaneous value of u, but rather only provide some mean values. That is, to 
employ (1) in practice usually requires an averaging or filter procedure. For this purpose, u is 
separated according to 

u = u + u , 

where u and u' refer to the averaged and the deviation quantities, respectively. The precise 
form of the averaging or filter method used determines additional calculation rules, e.g., ab = 
ab + a'b' for the classical Reynolds averaging. At the present stage it is not essential to already 
commit oneself to a definite averaging method. For nonlinear system (1) averaging usually gives 
expressions 

A'(j;,n(„),'u;) = 0, / = l,...,m, (2) 

where A' are smooth functions of their arguments whose explicit form is precisely determined 
by the form of A' and the chosen averaging rule. The tuple w = {w^,...,w^) includes all 
averaged nonlinear combinations of terms, which cannot be obtained by means of the quantities 
These combinations typically include such expressions as u'u', u'u, u'u'^, etc., referred to 
as subgrid scale terms. Stated in another way, system (2) contains more unknown quantities 
than equations. To solve system (2), suitable assumptions on w have to be made. An adequate 
choice for these assumptions is the problem of parameterization. 

The most straightforward way to tackle this issue is to directly express the unclosed terms w 
as functions of the variables x and u^^j.) some r which can be greater than n. In other words, 
system (2) is closed via 

A'(2;,n(„),/(2;,n(r))) = 0, l = l,...,m, (3) 

using the relation = f^{x, ■U(r)), s = 1, . . . , k. The purpose of this paper is to discuss different 
paradigms for the choice of the functions / = (/^, . . . , f^) within the symmetry approach, where 
k is the number of unclosed terms which are necessary to be parameterized. In other words, we 
should carry out, in different ways, group analysis of the class (3) with the arbitrary elements 
running through a set of differential functions. To simplify notation, we will omit bars over the 
dependent variables in systems where parameterization of w is already applied. 
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Remark. In the theory of group classification, any class of differential equations is considered 
in a jet space of a fixed order. That is, both the explicit part of the expression of the general 
equation from the class and the arbitrary elements can be assumed to depend on derivatives up 
to the same order. In contrast to this, for the construction of parameterization schemes it is 
beneficial to allow for varying the orders of arbitrary elements while the order of the explicitly 
resolved terms is fixed. This is why we preserve different notations for the orders of derivatives 
in the explicit part of the expression of the general equation and in the arbitrary elements of 
the class (3). 

2.1 Parameterization via inverse group classification 

Parameterizations based on Lie symmetries appear to have been first investigated for the Navier- 
Stokes equations. It was gradually realized that the consideration of symmetries plays a key role 
in the construction of subgrid scale models for the Navier-Stokes equations to allow for realistic 
simulations of flow evolution. See [29, 30] for a further discussions on this subject. The approach 
involving symmetries for the design of local closure schemes, was later extended in [38, 39, 40] 
in order to incorporate also the second law of thermodynamics into the consideration. 

For an arbitrary system of differential equations, this approach can be sketched as follows: 
First, determine the group of Lie symmetries of the model to be investigated. For common 
models of hydro-thermodynamics these computations were already carried out and results can 
be found in collections like [12]. Subsequently, determine the differential invariants of this group. 
If the left hand side of system (2) is formulated in terms of these invariants by an adequate choice 
of the function /, it is guaranteed that the parameterized system will admit the same group of 
point symmetries as the unfiltered system. Usually this leads to classes of differential equations 
rather than to a single model. That is, among all models constructed this way it is be possible 
to select those which also satisfy other desired physical and mathematical properties. 

The procedure outlined above can be viewed as a special application of techniques of in- 
verse group classification. Inverse group classification starts with a prescribed symmetry group 
and aims to determine the entire class of differential equations admitting the given group as 
a symmetry group [34]. Thus, in [29, 30, 38, 39, 40] it is assumed that the closure scheme 
for the subgrid scale terms leads to classes of differential equations admitting the complete Lie 
symmetry group of the Navier-Stokes equations. From the mathematical point of view, this 
assumption is justified as filtering (or averaging) of the Navier-Stokes equations introduces a 
turbulent friction term among the viscous friction term that already appears in the unfiltered 
equations. That is, filtering does not principally perturbs the structure of the Navier-Stokes 
equations. However, this assumption may not be as well justified if a model is chosen, where fil- 
tering leads to terms of forms not already included in the unfiltered model. In such cases, it may 
be more straightforward to solve the parameterization problem by inverse group classification 
only with respect to particular subgroups of the symmetry group of the initial model. 

The approach of inverse group classification usually relies on the notion of differential invari- 
ants [31, 34]. Differential invariants are defined as the invariants of the prolonged action of a 
given symmetry group. They can be determined either with the infinitesimal method [11, 34] 
or with the technique of moving frames [6, 9, 32]. In the present paper we will use the former 
method which is briefiy described here for this reason. 

Let X be the p-dimensional space of independent and U be the g-dimensional space of de- 
pendent variables. The connected Lie group G acts locally as a point transformation group 
on the space = X x U, with g denoting the associated Lie algebra of infinitesimal gen- 
erators. (The whole consideration is assumed local). Each element of q is of the form Q = 
^^{x,u)dxi + ip^{x,u)dua. In this section the indices i and j run from 1 to p while the indices a 
and b run from 1 to q, and the summation convention over repeated indices is used. The space 

= X X U(^r) is the rth prolongation of the space X x U (the rth order jet space), which is 
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the space endowed with coordinates Xi and u^, \a\ := ai + • • • + Op < r, where u'^ stands for 
the variable corresponding to the derivative d^°'^u°' /dx"^ ...dxp'', and a = {ai, . . . ,ap) is an 
arbitrary multiindex, Oj G N U {0}. The action of G can be extended to an action on and so 
the elements of g can be prolonged via 

Qir) = Q + E '^""^<' ■= or • • • Dp^'/'" - + e<+s.- (4) 

Here Dj = D^;- denotes the operator of total differentiation with respect to the variable Xj, i.e. 
Di = +Ua+5 '^u'^^ where 6i is the multiindex whose ith entry equals 1 and whose other entries 
are zero. More details can be found in the textbooks [31, 33, 34]. 

A differential function / (i.e., a smooth function from to M for some r) is called an (rth 
order) differential invariant of the group G if for any transformation g: {x,u) i— t- {x,u) from G 
we have that f{x,U(j.-^) = f(x,U(^r))- The function / is a differential invariant of G if and only if 
the equality Q(^r)f = holds for any Q E g. A vector field c) defined in J°° is called an operator of 
invariant differentiation for the group G if the result c)/ of its action to any differential invariant 
/ of G also is a differential invariant of G. 

The Fundamental Basis Theorem states that any finite-dimensional Lie group (or, more 
generally, any Lie pseudo-group satisfying certain condition) acting on possesses exactly 
p operators of invariant differentiation, which are independent up to linear combining with 
coefficients depending on differential invariants, and a finite basis of differential invariants, i.e., 
a finite set of differential invariants such that any differential invariant of the group can be 
obtained from basis invariants by a finite number of functional operations and actions by the 
chosen independent operators of invariant differentiation. 

For a vector field in J°° to be an operator of invariant differentiation of G, it is sufficient that 
it commutes with every infinitely prolonged operator from the corresponding Lie algebra g, i.e., 
[d, Q(oo)] = for any Q G g. If the group G is finite dimensional, a set of p independent operators 
of invariant differentiation can be found in the form d = /i*Dj by solving, with respect to the 
differential functions /i* = K'{x^U(^r))i the system of ffist-order quasi- linear partial differential 
equations 

where Q runs through a basis of the corresponding Lie algebra g and r equals the minimum 
order for which the rank of the prolonged basis operators of g coincides with its dimension. 
Eventually, it may be convenient to determine in the implicit form Q^{h^, . . . , K^, x, U(j.)) = 0, 
where det(r2;^,) ^ and $7-' satisfy the associated system of homogeneous equations 

(Q(,) + ih''Bi,e)dx^) = 0. 

In the infinite dimensional case, the construction of invariant differentiation operators is analo- 
gous though more sophisticated. 

A systematic approach to parameterization via inverse group classification hence consists of 
determining the basis differential invariants of a group together with the list of operators of 
invariant differentiation. Subsequently, there are infinitely many parameterizations that can be 
constructed, which admit the given group as a symmetry group. 

2.2 Parameterization via direct group classification 

The main assumption in the approach presented in [29, 40] is that a realistic subgrid scale model 
for the Navier-Stokes equations should admit the symmetry group of the original equations. 
However, this assumption is rather restrictive in more general situations. While it is true that a 
filtered model should be a realistic approximation of the unffitered equations, parameterization 
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schemes also have to take into account physical processes for which we may not have a precise 
understanding yet. That is, one eventually has to face the problem to deal with processes for 
which we may not even have a differential equation. This particularly means that a fixed set of 
symmetries (as for the Navier-Stokes equations) may not be obtainable. 

On the other hand, symmetries do provide a useful guiding principle for the selection of physi- 
cal models. As nature tends to prefer states with a high degree of symmetry, a general procedure 
for the derivation of symmetry- preserving parameterization schemes seems reasonable. The only 
crucial remark is, that we may not know in advance, which symmetries are most essential for 
capturing the characteristics of the underlying physical processes. For such problems, applica- 
tion of inverse group classification techniques is at once limited. Rather, it may be beneficial 
to derive parameterization schemes admitting different symmetry groups and subsequently test 
these various schemes to select among them those which best describe the processes under con- 
sideration. That is, instead of expressing the tuple w in system (2) using differential invariants 
of a symmetry group of the unfiltered equations (or another convenient symmetry group) from 
the beginning, we investigate symmetries of system (3) for different realizations of the functions 
/ which are eventually required to satisfy some prescribed conditions. This way, we could be 
interested in special classes of parameterizations, such as, e.g., time- or spatially independent 
ones. This naturally leads back to the usual problem of direct group classification: Let there 
be given a class of differential equations, parameterized by arbitrary functions. First determine 
the symmetries admitted for all choices of these functions, leading to the kernel of symmetry 
groups of the class under consideration. Subsequently, investigate for which special values of 
these parameter- functions there are extensions of the kernel group [34, 36]. 

To systematically carry out direct group classification, it is necessary to determine the equiva- 
lence group of the class, i.e. the group of transformations mapping an equation from the class (3) 
to an equation from the same class. Classification of extensions of the kernel group is then done 
up to equivalence imposed by the equivalence group of the class (3). 

The continuous part of the equivalence group can be found using infinitesimal methods in 
much the same way as Lie symmetries can be found using the infinitesimal invariance criterion. 
This firstly yields the equivalence algebra, the elements of which can then be integrated to give 
the continuous equivalence group. See [34, 36] for more details on this subject. 

We now formalize the method reviewed in the previous paragraphs. Let there be given a 
class of differential equations of the form (3), A'(x, /(x, M(r))) = Oi ^ = l,---,m. The 
arbitrary elements / usually satisfy an auxiliary system of equations S{x, U(^r)^ f{p){x, U(j.))) = 0, 
S = {S^, . . . , 5*), and an inequality S(x, f(p){x, ii(r))) = 0, where /(p) denotes the collection 
of / and all derivatives of / with respect to the variables x and U(^r) up to order p. The 
conditions S = and S 7^ restrict the generality of / and hence allow the design of specialized 
parameterizations. We denote the solution set of the auxiliary system by S, the system of 
form (3) corresponding to an / E 5 by Cf and the entire class of such system by C\s- 

The set of all point transformations that map a system Cj to a system Cj, where both 
/, / G 5 is denoted by T(/, /) and is referred to as the set of admissible transformations from 
the system Cf to the system Cj. The collection of all point transformations relating at least 
two systems from the class C\s gives rise to the set of admissible transformations of C\s- 

Definition 1. The set of admissible transformations of the class C\s is the set T{C\s) = 
{(/,/,^)l/,/G5,v5GT(/,/)}. 

That is, an admissible transformation is a triple, consisting of the initial system (with ar- 
bitrary elements /), the target system (with arbitrary elements /) and a mapping ip between 
these two systems. 

The usual equivalence group G~ = G^{C\s) of the class C\s is defined in a rigorous way in 
terms of admissible transformations. Namely, any element ^ from G~ is a point transformation 
in the space of (x,U(r),/), which is projectable on the space of {x,U(^r')) for ^^y < r' < r, so 
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that the projection is the r'th order prolongation of the projection of $ on the variables 

{x,u), and for any arbitrary elements / S 5 we have that $/ G S and S T(/, $/). The 

admissible transformations of the form if,^f,^\(x,u))^ where f £ S and ^ £ are called 
induced by transformations from the equivalence group G~ . Needless to say, that in general not 
all admissible transformations are induced by elements from the equivalence group. Different 
generalizations of the notion of usual equivalence groups exist in the literature [26, 36]. By 0~ 
we denote the algebra associated with the equivalence group G~ and call it the equivalence 
algebra of the class C\s- 

After clarifying the notion of admissible transformations and equivalence groups, we move on 
with the description of a common technique in the course of group analysis of differential equa- 
tions, namely the algebraic method. Within this method one at first should classify inequivalent 
subalgebras of the corresponding equivalence algebra and then solve the inverse group classifi- 
cation problem for each of the subalgebras obtained. This procedure usually yields most of the 
cases of extensions and therefore leads to preliminary group classification (see, e.g., [13, 14, 44] 
for applications of this technique to various classes of differential equations). 

The algebraic method rests on the following two propositions [7] : 

Proposition 1. Let a be a subalgebra of the equivalence algebra g~ of the class C\s, a C 
and let /'^(x, u^^)) £ S be a value of the tuple of arbitrary elements f for which the algebraic 
equation f = is invariant with respect to a. Then the differential equation C\fo is 

invariant with respect to the projection of a to the space of variables {x,u). 

Proposition 2. Let Si be the subset of S that consists of all arbitrary elements for which 
the corresponding algebraic equations are invariant with respect to the same subalgebra of the 
equivalence algebra g~ and let Oj be the maximal subalgebra of 0~ for which Si satisfies this 
property, i = 1,2. Then the subalgebras ai and 02 are equivalent with respect to the adjoint 
action of G~ if and only if the subsets Si and S2 are mapped to each other by transformations 
from G~. 

The result of preliminary group classification is a list of inequivalent (with respect to the 
equivalence group) members Cf oi the class £1^, admitting symmetry extension of the kernel of 
symmetry algebras using subalgebras of the equivalence algebra. 

Although the algebraic method is a straightforward tool to derive cases of symmetry exten- 
sions for classes of differential equations with arbitrary elements, there remains the important 
question when it gives complete group classification, i.e., preliminary and complete group clas- 
sification coincide. This question is of importance also for the problem of parameterization, as 
only complete group classification will lead to an exhaustive description of all possible parame- 
terization schemes feasible for some class of differential equations. The answer is that the class 
under consideration should be weakly normalized in infinitesimal sense, i.e., it should satisfy 
the following property: The span of maximal Lie invariance algebras of all equations from the 
class is contained in the projection of the corresponding equivalence algebra to the space of 
independent and dependent variables, 

(0/ 1 / G 5) c pr • 

At the same time, it is better to use a stronger notion of normalization introduced in [36]. 

Definition 2. The class of equations £1^ is normalized if V (/, /, ip) G T(Cs) 3<I> G G~: / = $/ 

and If = 

The normalization of £1^ in the sense of Definition 2 additionally implies that the group 
classification of equations from this class up to G~ -equivalence coincides with the group classifi- 
cation using the general point transformation equivalence. Due to this fact we have no additional 
equivalences between cases obtained under the classification up to G~ -equivalence. As a result. 
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solving the group classification problem for normalized classes of differential equations is espe- 
cially convenient and effective. 

In turn, depending on normalization properties of the given class (or their lacking), differ- 
ent strategies of group classification should be applied [36]. For a normalized class, the group 
classification problem is reduced, within the infinitesimal approach, to classification of subal- 
gebras of its equivalence algebra [3, 21, 36, 47]. A class that is not normalized can eventually 
be embedded into a normalized class which is not necessarily minimal among the normalized 
superclasses [36, 37]. One more way to treat a non-normalized class is to partition it into a 
family of normalized subclasses and to subsequently classify each subclass separately [36, 46]. If 
a partition into normalized subclasses is difficult to construct due to the complicated structure 
of the set of admissible transformations, conditional equivalence groups and additional equiv- 
alence transformations may be involved in the group classification [15, 35, 46]. In the case 
when the class is parameterized by constant arbitrary elements or arbitrary elements depending 
only on one or two arguments, one can apply the direct method of group classification based 
on compatibility analysis and integration of the determining equations for Lie symmetries up 
to G~ -equivalence [1, 28, 34]. Recall that these determining equations involve both coefficients 
of a Lie symmetry operator of a system Lf and the corresponding tuple of arbitrary elements / 
and follow from the infinitesimal invariance criterion [31, 34], 

Q(r')A'(x,'U(„),/(x,'U(r)))|£^. =0, / = 1, . . . ,m. 

Here r' = max{n,r}, the prolongation of Q is defined by (4) and the symbol means 
that above relation holds on solutions of the system £j. 



3 Symmetry preserving parameterizations 
for vorticity equation 

The inviscid barotropic vorticity equations in Cartesian coordinates reads 

Ct + W,C} = (5) 

where {a, 6} = axhy — aybx denotes the usual Poisson bracket with respect the variables x and y. 
The vorticity C and the stream function ip are related through the Laplacian, i.e. C = V^V- 
The two-dimensional wind field v = (u, v, 0)"'" is reconstructed from the stream function via the 
relation v = k x Vip, where k is the vertical unit vector. 

The maximal Lie invariance algebra go of the equation (5) is generated by the operators 

Vi =tdt-ipd^, dt, V2 = xdx + ydy + 2ipd^, 

J = -ydx + xdy, J* = -tydx + txdy + \{x^ + y'^)d^, 

X{-f^) = ^\t)dx - il{t)yd^, 3^(7') = l\t)dy + ^l{t)xd^, 

where 7"*^, 7^^ and x ^^un through the set of smooth functions of t. See, e.g. [2, 4] for further 
discussions. 

Reynolds averaging the above equation leads to 



+ ^,0 = v-(v'C0. (7) 

The term v'C = (ii'C, f'^', 0)^ is the horizontal eddy vorticity flux. Its divergence provides 
a source term for the averaged vorticity equation. The presence of this source term destroys 
several of the properties of (5), such as, e.g., possessing conservation laws. In this paper we aim 
to flnd parameterizations of this flux term, which admit certain symmetries. 
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A simple choice for a parameterization of the eddy vorticity flux is given by the down- gradient 
ansatz 

= -KVC, 

where the eddy viscosity coefficient K still needs to be specified. Physically, this ansatz accounts 
for the necessity of the vorticity flux to be directed down-scale, as enstrophy (integrated squared 
vorticity) is continuously dissipated at small scales. Moreover, this ansatz will lead to a uniform 
distribution of the mean vorticity field, provided there is no external forcing that counteracts this 
tendency [23]. The simplest form of the parameter K is apparently K = K{x,y), i.e. the eddy 
viscosity coefficient is only a function of space. More advanced ansatzes for K assume dependence 
on which is the eddy enstrophy [23] (see also the discussion in the recent paper [24]). 
This way, the strength of the eddy vorticity flux depends on the intensity of two-dimensional 
turbulence, which gives a more realistic model for the behavior of the fluid. There also exist a 
number of other parameterization schemes that can be applied to the vorticity equation, such as 
methods based on statistical mechanics [18] or the anticipated potential vorticity method [41, 45]. 

In the present framework, we exclusively focus on flrst order closure schemes. This is why 
we are only able to parameterize the eddy vorticity flux using the independent and dependent 
variables, as well as all derivatives of the dependent variables. This obviously excludes the 
more sophisticated and recent parameterization ansatzes of geophysical fluid dynamics from the 
present study. On the other hand, the basic method of invariant parameterization can already 
be demonstrated for this rather simple model. Indeed, symmetries of the vorticity equation 
employing the down-gradient ansatz or related parameterizations are investigated below using 
both inverse and direct group classiflcation. Physically more advanced examples for param- 
eterizations can be constructed following the methods outlined in Section 2 and exemplified 
subsequently. 

3.1 Parameterization via inverse group classification 

This is the technique by [29, 40] applied to the inviscid vorticity equation. In view of the descrip- 
tion of Section 2.1 this approach consists of singling out subgroups (subalgebras) of the maximal 
Lie invariance group (algebra) of the vorticity equation and computation of the associated differ- 
ential invariants (via a basis of differential invariants and operators of invariant differentiation) . 
These differential invariants can then be used to construct different parameterizations of the 
eddy vorticity flux. 

It is important to note that singling out subgroups of the maximal Lie invariance group of 
the vorticity equation is a meteorological way of group classification. This is why it is neces- 
sary to have a basic understanding of the processes to be parameterized before the selection 
of a particular group is done (otherwise, we would have to face the problem of how to com- 
bine these invariants to physically meaningful parameterizations). For the vorticity equation, 
we demonstrate the basic mechanisms of parameterizations via inverse group classification by 
singling out subgroups that allow to include the down-gradient ansatz of the previous section. 
This choice is of course not unique as there exist various other possibilities for parameterizations 
of the eddy vorticity flux. However, this choice allows us to demonstrate several of the issues of 
parameterization via inverse group classification. 

Invariance under the whole Lie symmetry group. To present differential invariants of the 
whole Lie invariance algebra Qq, we use the notation 

C — '4'xxx '^'4'xyyi ? — '^'4'xxy ^yyy^ ^ ~ ~^ V'xDy '4'y^x- 

The algebra go possesses no differential invariants up to order two. At the same time, it has 
the singular second order manifold determined by the equations ^ipxx = ipyy and ipxy = 0, which 
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is not essential for our consideration. A basis Bo of differential invariants of go consists of the 
third order differential invariants 

vc eve + T]Vr] {ve + 2rjC)^ + [Vr] - 2ecf 

Cx2 + Cy2' (^2 + ^2)1/2(^^2 + ^^2)' (^^2 + ^^2)3/2 

A complete set Oq of independent operators of invariant differentiation for this algebra is formed 
by the operators 

(02 ^ ^2)-l/2^^ (^^2 ^ ^^2)-l/2(^^j3^ ^ ^^D^)^ (^^2 ^ ^^2)-l/2 j^^ _ ^^j^^y 

The computation of Bq and Oq is cumbersome and will be presented elsewhere. At the same 
time, the result of the computation can be checked in a rather direct and simple way. Indeed, 
the cardinality of Oq equals three. The elements of Oq are linearly independent over the ring 
of differential invariants of go and commute with the infinite prolongations of all vector fields 
from the generating set (6) of go- Since each element I of Bq satisfies the condition (5(3)/ = 0, 
where the operator Q runs through the operators (6), it is a differential invariant of go. The 
invariants belonging to Bq are functionally independent. Moreover, for any fixed order r an rth 
order universal basis of differential invariants of go can be constructed via acting by operators 
from Oq on invariants from Bq. We only sketch the proof of the last assertion. The cardinality 
of any rth order universal basis of differential invariants of go, where r ^ 4, equals the difference 
between the dimension of the jet space and the rank of the rth prolongation of go, 

A^ = 3+(' + ^)-(3r + 8). 

Acting on elements of Bq by operators from Oq k — 'i times, 3 ^ A; ^ r, we obtain a set of 
/c-order differential invariants which is of maximal rank with respect to the /c-order derivatives 
involving at least two differentiations with respect to space variables. Choosing, for each fe, 
a subset of invariants associated with a nonzero /c-order minor in the corresponding Jacobi 
matrix and uniting such subsets for k ^ r, we construct exactly N functionally independent 
differential invariants of order not greater than r, which hence form an rth order universal basis 
of differential invariants of go. 

The above case where invariance of the parameterization under the whole symmetry group 
of the vorticity equation is desired can be neglected for physical reasons. This is since it is 
impossible to realized, e.g. the down-gradient ansatz described in the previous section within 
this framework. It can easily be checked that the corresponding vorticity equation with parame- 
terized eddy vorticity fiux only admits one scaling operator for any physically meaningful ansatz 
for K. In contrast to the example of the Navier-Stokes equations discussed in [29], the vorticity 
equation hence does not allow physical parameterizations leading to a closed model invariant 
under the same symmetry group as the original vorticity equation. This is why it is beneficial 
to single out several subgroups of the maximal Lie invariance group and consider the invariant 
parameterization problem only with respect to these subgroups. 

Explicit spatial dependency. If the two-dimensional fluid is anisotropic and inhomogeneous 
the only subalgebra of (6) that can be admitted is spanned by the operators 

For this subalgebra, a basis of invariants is formed by x, y, ipx and ipy. Independent operators 
of invariant differentiation are exhausted by Dj, D^. and D^^. If we express the right hand side 
of (7) in terms of differential invariants of the above subalgebra, a possible representation reads 

Ct + {i^,C} = K{x,y)V'C. 
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Hence we assembled our parameterization using the (differential) invariants x, y, T)^ipx = '4^xxxx, 
DyDxTpx = "ipxxyy and ^yipy = ''Pyyyy This boils down to the usual gradient ansatz for the eddy 
flux term, where the eddy viscosity K explicitly depends on the position in the space. Note, 
however, that this ansatz is only one possibility which is feasible within this class of models. 

Rotationally invariant fluid. In case the two-dimensional fluid is isotropic, the resulting 
parameterized system should also admit rotations. Hence, we seek for differential invariants of 

the subalgebra j spanned by the operators 

dt, J = xdy-ydx, = txdy~tydx + \{x^ + y^)d^,, = xW^V- 

A complete set of independent operators of invariant differentiation is formed by V , xTi^ + J/Dy 
and — yDx + xDy and a basis of differential invariants consists of 

r = y/x"^ + xiljy - yVx, r'^ii^xx + ijjyy) - '2{xijJx + vi^y), 

r^{xiptx + yipty) + {xipx + yil)y){xy{'il}yy - ipxx) + {x^ - y^)'4^xy + xil^y - yipx)- 

In the polar coordinates (r, (^) with (p = arctany/x, these sets have, after an additional rear- 
rangement, simpler representations O = {Df, D^, D<^} and B = {/", a = 0, . . . , 3}, respectively, 
where Dj = Dt -|- r''^ipr^ip and 

Any element of O indeed is an invariant differentiation operator for j since it commutes 
with the infinite prolongation of every vector field from j. The fact that 1°', a = 0, . . . , 3, are 
differential invariants of j is also checked in a rather direct way, by the substitution to the 
condition Q(2)I = 0, where the operator Q runs through j. These invariants obviously are 
functionally independent. 

The most difficult part is to prove that for any fixed order r we can construct an rth order 
universal basis of differential invariants of j by invariant differentiations of 7°, a = 0, ... ,3. The 
number of elements in any rth order universal basis of differential invariants of j, where r ^ 1, 
equals 

3+(^^^)-(^ + 4) = ^(r + l)(r + 5) 

(the dimension of the jet space minus the rank of the rth prolongation of j). The commutation 
relations between the operators of invariant differentiation are 

[D^, D^] = 0, [D^, Dt] = (r- Vr)rD^, [D^, Dt] = r" Vr^.D^'- 

The elements 7", a = 1, . . . , 3 of H can be represented in the form 7^ = D^^, 7^ = Dj.r~^Dj.^ 
and 7^ = DtD^V'- Hence, acting by the operators of invariant differentiation on elements of B, 

we can construct 

, /r-H-3\ fr-2 + 2\ /r - 2 + l\ r, 

functionally independent invariants of order not greater than r (the zeroth order invariant r plus 
acting on 7^ by any operators from B at most r — 1 times plus acting on 7^ by Dt or D^ at most 
r — 2 times and plus acting on 7^ by D^ at most r — 2 times). In view of coinciding the above 
numbers, the proof is completed. 

As r, + {V') C} ^'^d V^C are differential invariants of j, the parameterizations of the form 

are invariant with respect to j. 

In the same fashion it would be possible to derive classes of parameterizations that preserve 
other subalgebras of Qq, e.g. including (generalized) Galilean symmetry or a scaling symmetry, 
but we do not derive them in this paper. 
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3.2 Equivalence algebras of classes of generalized vorticity equations 

In order to demonstrate different possible tecliniques, we present tlie details of the calculation 
of the usual equivalence algebra Qi for the class of equations 

Ct + {i^,C} = ^if{t,x,y,C.Xy)=fi+flQj, C:=^«, (8) 

where for convenience we introduce another notation for the independent variables, t = zq, 
X = zi and y = Z2, and omit bars over the dependent variables. Throughout the section the 
indices i, j and k range from 1 to 2, while the indices k, A, fi and v run from to 2. The 
summation over repeated indices is understood. A numerical subscript of a function denotes the 
differentiation with respect to the corresponding variable z^. 

In fact, the equivalence algebra of class (8) can be easily obtained from the much more general 
results on admissible transformations, presented in Section 3.3. At the same time, calculations 
using the direct method applied for finding admissible transformations are too complicated and 
lead to solving nonlinear overdetermined systems of partial differential equations. This is why 
the infinitesimal approach is wider applied and realized within symbolic calculation systems. 
The usage of the infinitesimal approach for the construction of the equivalence algebra of (8) 
has specific features richly deserving to be demonstrated here. 

Theorem 1. The equivalence algebra of class (8) is generated by the operators 

Vi = tdt - i^d^ - Cx^c, - Cydc, - 2fdfi - 2fdp , dt, 

T>2 = xd, + ydy + 2i;d^ - G% - Cy^c, + fOp + fdp , 

J{P) = Pxdy - ^yd, + f + y^)d^ + /5(Cx5c, - CydO 

+ {Pttx - f3f)dji + i(3uy + Pf)df2, (9) 
^(V ) = I'd, - j^yd^, 3>(7') = l^dy + ^fxd^, 
■^(f^) = ^{x^ + y^){d^p + Cydfi - Cxdp) + atxdfi + atydp, 
i{(S) = 6{d^ + Cydfi - Cxdf2), g{p) = p^dp - pydfi, Z{x) = xd^p, 

where f3, 7*, a and x <ire arbitrary smooth functions oft solely, 5 = 6{t,x,y) is an arbitrary 
solution of the Laplace equation 6xx + Syy = and p = p{t, x, y) is an arbitrary smooth function 
of its arguments. 

Remark. Although the coefficients of and d(^y can be obtained by standard prolongation 
from the coefficients associated with the equation variables, it is necessary to include the corre- 
sponding terms in the representation of the basis elements (9) in order to guarantee that they 
commute in a proper way. 

Remark. The operators G{p) and Hix) — ^ix) arise due to the total divergence expression 
of the right hand side of the ffist equation in (8), leading to the gauge freedom in rewriting 
the right hand side of the class (8). They do not generate transformations of the independent 
and dependent variables and hence form the gauge equivalence subalgebra of the equivalence 
algebra (9) [36]. The parameter- function p is defined up to summand depending on t. 

Proof. As coordinates in the underlying fourth-order jet space ^^^^ , we choose the variables 

Zf,, V, V'm, ipf„y, ps^u, ^/^A^^, A ^ /i ^ z^, (/i, I/) / (2,2), Cm> 
K ^ A ^ /i ^ I/, {p,v)^ (2, 2), Q^y, pi^u. 

(Variables of the jet space and related values are defined by their notation up to permutation of 
indices.) The variable Co of the jet space is assumed principal, i.e., it is expressed via the other 
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coordinate variables (called the parametric ones) in view of equation (8). Under calculation we 
also carry out the substitutions ip22ii = C^l~'^ll^l■ To avoid repetition of the above conditions for 
indices, in what follows we assume that the index tuples (A, /i, i^) and (k, A,//, z^) satisfy 

these conditions by default. 

Due to the special form of the arbitrary elements we have to augment equation (8) with 
the following auxiliary system for /*: 

P = P = P = P = P = P = P = (^()) 

As we compute the usual equivalence algebra rather than the generalized one [26] and the 
arbitrary elements /* do not depend on fourth order derivatives of ip, the elements of the algebra 
are assumed to be vector fields in the joint space of the variables of 5*^^-' and the arbitrary elements 
/*, which are projectable to both the spaces {t,x,y,tp) and In other words, the algebra 

consists of vector fields of the general form 

Q = eo, + vd^ + v^d^^ + v^'d^^^ + v^^^^d^,^^ + e^d^^ + ^'df., 

where .^'^ = ^'^{t, x,y,^), rj = r]{t, x,y,^), the coefficients corresponding to derivatives of ip are 
obtained by the standard prolongation (4) from and r], the coefficients 9'^ are obtained by the 
standard prolongation from and 9 = rj^"^, and the coefficients 99* depends on all the variables 
of 2^^^ and the arbitrary elements f^. As a result, each element from the equivalence algebra 
is determined by its coefficients rj and ip'''. To act on the equations (8) and (10) by the 
operator Q, we should additionally prolong it to the variables ^JuXf^u and C^lu in the conventional 
way and to the derivatives of /, assuming all the variables of ^^^^ as usual ones: 

Q = Q + r^^^^^d^^,^^ + 9^'^dc^^^ 

First we consider the infinitesimal invariance conditions associated with equations (10). The 
invariance condition for the equation = is 

Splitting with respect to derivatives of /* in the latter equation implies that c/?^ = 0, = 0, 
9lp = 0. As 6** = DjDjD.j(7/ — C'^tp^) + ^^'^iijji, we additionally derive the simple determining 
equation rj^^ = 0. 

In a similar way, the invariance conditions for the equations /^^ = 0, /^^^ = 0, /^^ = 
and /^^ = can be presented in the form 

(io) = '^W-^'.4=0' 

^ lEq. (10) - '^^Am-. ~ '^^Xf.JCk - ^' 

^'■'°Ie,. {io)=V'Co-^Co4=0' 

which is split into v^;^ = 0, = 0; ^^^^^ = 0, 9^^ = 0; ^^^^^ = 0, 9^^^^ = 0; and = 0, 
9^^ = 0, respectively. The equations 9^^ = 0, 9^^^^ = and 9^^ = provide no essential 
restrictions on the coefficients S^'^, rj and From the equation 9^^ = we derive that = 0, 
^1 + ^2 ^ and Cl - = 0. Hence 

9 = n" = rjjj + 2??^^^^ + wV'ii - '^Cj'^ij = Vjj + '^Vj^P'^j + (w ~ '^^i)C- 
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It remains to solve the determining equations following from the invariance condition for 
equation (8). The invariance condition reads 

0" + v'C2 + V'l^' - r/'Ci + ^20' = + ip'^^Cji + fi/\ 

or explicitly 

HV^ - 2^1 - + flC^J - i^lC2 + i^2Cl) - iK^ 

+ivi + Vi(>'^i - ?iV'i)C2 + ipii'njj2 + Vjji^'^2 + 2r?2j^V'j + 2??jvV'2i + ~ 2^1 )C2 - ^20) 
-{m + r]^ip2 - C2A)Ci - 4^2{r]jji + ^jji^i^i + 2??ij^V'i + 2??jv^ii + iVi' - 2d)Ci - CiCi) 

= ^ + ^pfl - siif] - o'Jl + OM^ + V'^iJ - olfj + flfK 

Collecting the coefficients of gives rjj^ = 0. This implies that 6*^ = 0. Similarly, the 
coefficients of ipiQj lead to the equation rnjj = and r/^ — + = 0. As = and rn^ = 0, 
the second equation together with the relations ^1 = ^1 ^2 + = implies that = 0. 
Then, the coefficient of C, gives = and the coefficients of /j lead to (^^2 = ^2) V'/i = ^''^d 
(/j^i = (^^2 = ~ 2^^. In view of the determining equations that we have already derived, the 
terms involving are identically canceled. Note that the coefficients of f^ Cki simultaneously 
lead to the same set of equations as the coefficients of /j. 

The remaining part of the invariance condition is 7]jjt — + '?iC2 — ^2^1 = ^\ + dj^^j- 
Splitting with respect to dj in this relation gives ip^^ = ip'^^ = 0, ip^^ + ip'^^ = and 

fi = Vjjt - CtCi + mC2 - mCi- 

Acting on the last equation by the operator djdi^., we obtain = 0. Further splitting with 
respect to Ci and C2 is not possible since ip^ may depend on them. 
Finally, the reduced system of determining equations reads 

e°=^° = e^i = o, e; = ^fc = o, 4 = 0, ei = f2, ^2+^1 = 0, 

= 0' Vi^ = 2^1 - Vijj = 0, 
^;> = 0' ^W=0' ^W.=0' '/'k..=0' ^Co=0' 

,^1 fl ,o2 t2 ,^1 ,^2 tl otO 

¥'/2=?2' 9^/1= ?ii V'/i = = ?i - , 

The solution of this system provides the principal coefficients of the operators from the equiva- 
lence algebra of the class (8): 

f = cit + co, ^1 = C2X - /3y + 7\ f = Px + C2y + , 

(J 

77 = (2c2 - ci)V + d- ily + ijx + ^{x^ + y^) + -{x^ + y^) + x, 

a (11) 
= (C2 - 2ci)/i - 13 f + 6Cy + 2 (^^ + y^)Cs/ + + - py, 

= Pf + (C2 - 2ci)/2 - 6C. - f + y')Cx + Aty + aty + p., 

where /3, 7*, o" and x are real- valued smooth functions of t only, cq, ci and C2 are arbitrary 
constants, p is an arbitrary function of t, x and y and 6 = 6{t, x, y) is an arbitrary solution of 
the Laplace equation 5jj = 0. 

Splitting with respect to parametric values in (11), we obtain the coefficients of the basis 
operators (9) of the algebra . Recall that the coefficients 77^, i]^^ , 7^-**/'^ and 9'^ are calculated 
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from and 7] via the standard procedure of prolongation and the coefficients do not depend 
on Tpfj^, ip^,^, tpx^u, and Co- Therefore, both the operators from g^" and their commutators are 
completely determined by the coefficients of d^, d^, 9^. and dp. This is why in (9) and similar 
formulas we omit the other terms for sake of brevity. □ 

Remark. The auxiliary system for the arbitrary elements is an important component of the 
definition of a class of differential equations. Its choice is usually guided by some prior knowl- 
edge about the processes to be parameterized. We have decided to assume that the arbitrary 
elements and depend also on t, keeping in mind two more, purely mathematical, rea- 
sons. The first reason is that the projection of the corresponding equivalence algebra on the 
space (t, X, y, ip) contains the maximal Lie invariance algebra go of the vorticity equation (5) 
which is the initial point of the entire consideration. The basis operators (6) of go are obtained 
from (9) by 

Pi=PPi, dt = Vdt, V2 = PV2, J = PJ{1), j' = PJit), 

Af(y) = p^(7i), y{j^) = Fyij^), zix) = P^(x), 

where P denotes the projection operator on the space {t,x,y,Tp). (Though the expressions 
for the operator dt (resp. Af(7^), 3^(7^) or Z{x)) and its projection formally coincide, they in 
fact determine vector fields on different spaces.) The second reason is that the class (8) is 
normalized, cf. Section 3.3. This in particular implies that the maximal Lie invariance algebra 
of any equation from the class (8) is contained in the projection of the equivalence algebra g^" 
of this class. 

We also calculate the equivalence algebras of two subclasses of the class (8). 
The first subclass corresponds to parameterizations not depending on time explicitly and, 
therefore, is singled out from the class (8) by the further auxiliary equation 

fi = 0, 

which has no influence on splitting of the invariance conditions for the equations (8) and (10) 
and gives the additional determining equations (pi = Q = 61 = 0. These determining equations 
imply that /3, 7* and a are constant, 5 is a function only of x and y and p can be assumed as a 
function only of x and y. Therefore, the equivalence algebra of this subclass is 

(Pi, dt, V2, J{i), y{i), n{i), n{5), g{p), z{x)), 

where the parameter-function 6 = 5{x, y) runs through the set of solutions of the Laplace 
equation 5xx + ^yy = and p = p{x, y) is an arbitrary function of its arguments. 

The second subclass is associated with spatially independent parameterizations. Hence we 
additionally set 

/; = o. 

It has to be noted that after attaching this condition we cannot split with respect to /j as we 
did in the course of solving the determining equations. However, precisely the same conditions 
obtained from splitting with respect to /j can also be obtained from splitting with respect to 
/^.. Hence the condition /j = only leads to the additional restriction ip^j = and, therefore, 
we find that (5j = 0, o" = 0, = and pij = 0. Without loss of generality we can set p = p'^{t)zi, 
where are arbitrary smooth functions of t. As a result, the equivalence algebra g^ of the 
second subclass is generated by the operators 

A, dt, V2, J{1), J{t), 3>(7'), n6), Gip'x + p^y), Z(x), 
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where 7', 5 and x arbitrary smooth functions of t. 

The intersection of the above subclasses corresponds to the set of parameterizations indepen- 
dent of both t and (x, y) and is singled out from the class (8) by the joint auxiliary system 

fi = f; = o. 

Its equivalence algebra is the intersection of the equivalence algebras of the above subclasses 
and, therefore, equals 

(Pi, dt, f>2, J{1), 3>(1), n{l), Gip'x + p^y), Z{x)), 

where p^, p^ and x ^'^^ arbitrary smooth functions of t. 

3.3 Normalized classes of generalized vorticity equations 

In the course of computing the set of admissible transformations of a class of differential 
equations, it is often convenient to construct a hierarchy of normalized superclasses for this 
class [36, 37]. This is why here we also start with the quite general class of differential equations 

Qt- F{t,x,y,'4),'4jx,'4^y,C,,CxXyXxx,C,xy,C,yy) C:=ipii, (12) 

where {F^^ , F,^^ , F^^^ , F^^^ , F/^^^ ) / (0, 0, 0, 0, 0) , to assure that the generalized vorticity equations 
of the form (8) belong to this class. We use notations and agreements from the previous section. 
In particular, z = {zq,zi,Z2) = {t,x,y), the indices i, j and k again run through {1,2}, while 
the indices k, A, p and range from to 2. 

Admissible transformations are determined using the direct method in terms of finite trans- 
formations. Namely, we aim to exhaustively describe point transformations of the form 

T: z^ = Z^(z,V), i^ = ^{z,i^), where J = ^^f ^ 0, 

d[Zo,Zi,Z2,'ll^) 

which map an equation from class (12) to an equation from the same class. We express deriva- 
tives of the "old" dependent variable V with respect to the "old" independent variables z via 
derivatives of the "new" dependent variable ip with respect to the "new" independent variables z. 
The latter derivatives will be marked by tilde over tp. Thus, the derivative of ip with respect 
to z^ is briefly denoted by tpfj,, etc. Then we substitute the expressions for derivatives into the 
equation Q — F = 0, exclude the new principal derivative V'022 using the transformed equation 
■0022 = —"0011 + F, split with respect to parametric variables whenever this is possible and solve 
the obtained determining equations for and ^ supplemented with the inequality J 7^ 0, 
considering all arising cases for values of the arbitrary element F and simultaneously finding the 
expression for F via F, and ^. 

The first order derivatives ip^ are expressed in the following manner: 

_ ^,-i..z; _ 

where we have introduced the notation V = V{z,jp,z) := ^(z,-!/;) — ipu{z)Z^ {z^ip) which is 
assumed as a function of the old dependent and independent variables and the new independent 
variables, so that V"^ = ^ ^ — ipyZ'l^ and = — ^yZ^. We will not try to express the old 
variables via the new variables by inverting the transformation. This is a conventional trick 
within the direct method, which essentially simplifies the whole consideration. In what follows 
we will also use three more abbreviations similar to V^: 

^i^v ._ ^My^ _ z^A^^, W^"" := U^^'U"^ Fq^ , pi" ■= U^"^ - V'Ft^. . 
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Higher order derivatives are expressible in an analogous way. The Laplacian of ■0i e.g., reads 



For the class (12) considered here, we need the derivatives of the Laplacian up to second order. 
The highest derivatives required are of the form 

i^iijk = V^'U^^'W^'U^^U^'i^^.^x + . . . , 

where the tail contains only derivatives of tp up to order three. 

Denote by G the left hand side of the equation obtained by substituting all the expressions for 
derivatives into (12). For the transformation T to be admissible, the condition G j = has 
to be satisfied for any tuple of the subscripts (fi, v, k, A) in which at least one of the subscripts 
equals 0. Under varying the subscripts, this condition leads to the following system: 



G - 

■00000 


= 0: 


G - 

V'OOOi 


= 0: 


VOOij 


= 0: 



00 



Suppose that [/O'^f/O^ / 0. Then the above equations imply that := U^'V^Fq. = 0. If 

rank(C/'^*) < 2 then for any ^ and v 



V (9(2;i,Z2,V') d{zi,Z2,'4>) J 



and after splitting with respect to tp\ we obtain that 



d{zi,Z2,i>) d{zi,Z2,'lp) 



or Z^ = ^^ = 0, 



but this contradicts the transformation nondegeneracy condition J ^ 0. Hence rank(C/'^') = 2 
and, therefore, the equation U^^U'^^ Fq. = sequentially implies that U^^ Fq. = and Fq. = 0. 
Then, the necessary conditions G^^^^ = and G^^^_ = for admissible transformations are 
respectively equivalent to the equations U^^U^^P^ = and U^^U^^P' + 2U^^U'^P^ = which 
jointly gives in view of the condition jjOkjjOk _^ q ^j^g^^ pM — q_ Thus, we should have det{U^'^) = 
0. At the same time, 

detiU'^n = V/{\Zl^, Zlzl\V^ - IK, Zl Zl\Zl - \Zl K, Zl\Zl - \Zl Zl V,\Zl) 

d{Z0,Zi,Z2,1p) 

that leads to a contradiction. Therefore, the supposition JJ^'^JJ^^ 7^ is not true, i.e., JJ^^U^^ = 
and hence U^^ = 0. Substituting the expressions for U^^ and V into the last equation and 
splitting with respect to V'ai, we derive the equations 

zlz^^ = zlzl z,% = 4m/,. 

The tuples (Z{',^'i), (Z2,^'2) and (Z^,^'^) are not proportional since J 7^ 0. This is why we 
finally obtain the first subset of determining equations Z^ = Z^ = 0. It follows from them that 
Zq ^ (otherwise J = 0) and expressions for "old" derivatives with respect to only x and y 
contain "new" derivatives only of the same type. In other words, derivatives of ip involving 
differentiation with respect to t appear only in the expressions for tpoaa and we can simply split 
with respect to them via collecting their coefficients. 
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Equating the coefficients of V'012 leads, in view of the condition Zq 7^ 0, to the equation 
jjikjj2k ^ 0, i.e., 

(4^1/^ - Zl^k + (4^1 - ZlZ^^)i>2) {Zl^^ - Zl^k - {Zlzl - ZlzDi^i) = 0. (13) 

We spht equation (13) with respect to and ^'2- Collecting the coefficients of ■0iV'2 gives 
the equation {Z^Zl - ZlzDiZ^Zl - ZlZ^ = 0, or equivalently Z^Z^ - ZlZ^ = 0. As 
rank(ZJ, Z2, Z^) = 2, this implies that = and, therefore, ^'^ / 0. Consequently, equa- 
tion (13) is reduced to Z^Z"^ = 0. 

The derivative "0022 is assumed principal, 1/^022 = — V'oii + F- Hence another third order 
derivative of the above type appropriate for splitting is only V'oii- The corresponding equation 
Z^Z^ = Z'^Z'^ := L joint with the equation Z^Z'^ = implies that the functions Z^ and Z'^ 
satisfy the Cauchy-Riemann system = eZ|, Z2 = sZf, where e = ±1, and hence = 0. 
Note that L / since J / 0. 

Analogously, collecting the coefficients of V'oi and further splitting with respect to ipj lead to 
the equations Z^Z-j^"^^^ = and Z^^^kip = 0. Therefore, = and "^k^ = 0. Here we take 
into account the inequalities L 7^ and det(Z^) / 0. 

We do not have more possibilities for splitting. The derived system of determining equations 
consists of the equations 

Zk = Z^ = 0, = 0, zlzl = 0, zlzl = zlzl, ^'^^ = "^ktp = 0. 

The remaining terms determine the transformation rule for the arbitrary element F. This is 
why any point transformation satisfying the above determining equations maps every equation 
from class (12) to an equation from the same class and, therefore, belongs to the equivalence 
group Gi of class (12). In other words, any admissible point transformation of class (12) is 
induced by a transformation from , i.e., class (12) is normalized. As a result, we have the 
following theorem. 

Theorem 2. Class (12) is normalized. Its equivalence group G^ consists of the transformations 
i=T{t), x = Z\t,x,y), y = Z\t,x,y), = T{t)tl; + ^{t,x,y), 

where T, Z^ , T and ^ are arbitrary smooth functions of their arguments, satisfying the conditions 
zlzl = 0, zlzl = zlzl := L, TtTL 7^ 0, and the subscripts 1 and 2 denote differentiation 
with respect to x and y, respectively. 

The expression for the transformed vorticity is also simple: Q — L ^(T(^ + $m). 

Remark. The continuous component of unity of the group Gi consists of the transformations 
from Gi with Tt > 0, e = 1 and T > 0. Therefore, a complete set of independent discrete 
transformations in Gi is exhausted by the uncoupled changes of the signs of t, y and ip. In 
particular, the value e = — 1 corresponds to alternating the sign of y. 

Consider the subclass of class (12), singled out by the constraints = 0, F^^ = —C,y and 
P-tpy = Cx) i-e., the class consisting of the equations of the form 

C,t + ^xCy-'<l^yC,x = H{t,X,y,C,C,x,C,y,C,xxXxy,Cyy), C ■= '<Pii, (14) 

where H is an arbitrary smooth function of its arguments, which is assumed as an arbitrary 
element instead oi F = H — tpxCy + ipyCx- The class (14) is still a superclass of the class (8). 
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Theorem 3. Class (14) is normalized. The equivalence group of this class is formed by the 
transformations 



t = T, X = X{xc — ys) + 7 , sy = X{x5 + yc) + 7 



2 / cr 



where e = ±1, c = cos/3, s = sin/3; r, A, /3, 7* and a are arbitrary smooth functions of t 
satisfying the conditions A > and Tt 7^ 0; 5 = 5{t, x, y) runs through the set of solutions of the 
Laplace equation bxx + ^yy = 0. 

Proof. The class (14) is a subclass of the class (12) and the class (12) is normalized. Therefore, 
any admissible transformation of the class (14) is generated by a transformation from the equiv- 
alence group of the superclass. It is only necessary to derive the additional restrictions on 
transformation parameters caused by narrowing the class. 

The group G"^ is a usual equivalence group [34], i.e., in contrast to different generalizations of 
equivalence groups [25, 36], it consists of point transformations of the joint space of the equation 
variables and arbitrary elements, and the components of transformations for the variables do not 
depend on the arbitrary elements. Any transformation from G^" is additionally projectable to the 
space of the independent variables and the space of the single variable t. This is why it already 
becomes convenient, in contrast to the proof of Theorem 2, to express the new derivatives via 
old ones. Then we substitute the expressions for new derivatives into the transformed equation 
Ct + ^iCy ~ V'yCi = exclude the principal derivative '4)tyy using the equation 

split with respect to parametric variables whenever this is possible and solve the obtained deter- 
mining equations. As equations from the class (14) involve derivatives ■0x and ipy in an explicitly 
defined (linear) manner, we can split with respect to these derivatives, simply collecting their 
coefficients. Since these coefficients do not involve the arbitrary element H, we can further split 
them with respect to other derivatives. As a result, we obtain the equations 

T = e^, Li = 0, ^^jj^ = 0, 

where e = ±1 and other notations are defined in the proof of Theorem 2. Therefore, L and 
are functions of t only. As L > 0, we can introduce the function A = ^/L of t. Acting by 
the Laplace operator djj on the conditions Z\Z\ = A^ and Z'^Z'^ = A^ and taking into account 
that Z' are solutions of the Laplace equation, = 0, we derive the important differential 
consequences Zj^ = 0, which imply that the functions are affine in (x,y). Hence there exists 
a function /3 = /3(t) such that Z\ = Ac and Zg = —As, where c = cos/3 and 5 = sin/3, and, 
therefore, Z\ = eAs and Z2 = eXc. We re-denote T by r for the sake of notation consistency 
and represent $ in the following form^: 

$ = 5{t, X, y) + ^(x^ + /) + e- (^^x'^ + y^) - ^l{x5 + yc) + ^"fixc - ys) 
2 Tt \2 



^There is an ambiguity in representations of and For example, the last summand in the representation 
of $ can be omitted. The usage of the above complicated representations is motivated by a few reasons: the 
consistency with the notation of basis operators of the equivalence algebra gi from Theorem 1, the simplification of 
the expression for the transformed arbitrary element H and the convenience of studying admissible transformations 
within subclasses of the class (14). 
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where o" is a function of t and 5 = 5{t, x, y) is a solution of the Laplace equation 6xx + 5yy = 0. 
Collecting the terms without and ipy gives the transformation for the arbitrary element H. 

Similarly to the proof of Theorem 2, any transformation from satisfying the above addi- 
tional constraints maps every equation from the class (14) to an equation from the same class 
and, therefore, belongs to the equivalence group of the class (14). In other words, any 
admissible point transformation of the class (14) is induced by a transformation from G^, i.e., 
the class (14) is normalized. □ 

Remark. The transformations from the equivalence group G^, which are associated with the 
parameter- function 6 depending only on t, and only such transformations identically act on the 
arbitrary element H and, therefore, their projections to the space of independent and dependent 
variables form the kernel (intersection) of point symmetry groups of the class (14). 

Corollary 1. The subclass of the class (14) singled out by the constraint H,^ = is normalized. 
Its equivalence group G3 consists of the elements o/G^ with tu = 0. 

Proof. As the vorticity and its derivatives are transformed by elements of G^ according to the 
formulas 



C = :^(C + /3t) + 2^, = :^0, (16) 

it follows from (15) under the constraints = and = that tu = 0. The rest of the proof 
is similar to the end of the proof of Theorem 3. □ 

Corollary 2. The subclass of the class (14) singled out by the constraints Hi = is normalized. 
Its equivalence group G4 consists of the elements o/G^ with At = 0, o" = and 6ij = 0. 

Proof. As any admissible transformation of the class (14) has the form (15) and, therefore, the 
vorticity and its derivatives are transformed according to (16), the system Hx = 0, Hy = is 
equivalent to the system Hx = 0, Hy = 0. After differentiating the last equation in (15) with 
respect to x and y and splitting with respect to (x and (y, we derive all the above additional 
constraints on transformation parameters. The rest of the proof is similar to the end of the 
proof of Theorem 3. □ 

Corollary 3. The subclass of the class (14) singled out by the constraints H(^ = and Hi = 
is normalized. Its equivalence group G^ consists of the elements of G2 with tu = 0, Xt = 0, 
cr = and 5ij = 0. 

Proof. The subclass under consideration is normalized as it is the intersection of the normalized 
subclasses from Corollaries 1 and 2. Therefore, we also have G^ = G3 H G4 ■ □ 

Remark. For the subclass from Corollary 3, the kernel of point symmetry groups is essen- 
tially extended in comparison with the whole class (14). It is formed by the projections of 
elements of the equivalence group G2 , associated with the parameter- functions 7^ and 7^ and 
the parameter-function 5 depending only on t, to the space of independent and dependent vari- 
ables, cf. Section 3.4. 

A further narrowing is given by the condition that the arbitrary element H with H(^ = is a 
total divergence with respect to the space variables, i.e., H = Dj/' for some differential functions 
/* = /*(^' ^) Cx) Cy)- The corresponding subclass rewritten in the terms of /* coincides with 
the class (8) and is singled out from the class (14) by the constraints H^ = and EH = 0, where 
E = — ^id^ + X^j^j + . . . is the associated Euler operator. In this Euler operator, 

the role of independent and dependent variables is played by {x, y) and C, respectively, and the 
variable t is assumed as a parameter. The vorticity C, can be considered in E as the dependent 
variable instead of il) since the arbitrary element H depends only on combinations of derivatives 
of V' being derivatives of Q. 
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Remark. It is obvious that the arbitrary element H satisfies the constraints = and EH = 
if it is represented in the form H = Dj/* for some differential functions /* = P{t, x, y, Qx^ Cy)- The 
converse claim should be proved. Thus, the constraint EH = implies the representation H = 
Dj/* for some differential functions f^{t, x,yX, Cx, Cy)^ which may depend on Substituting 
this representation into the constraint H(^ = and splitting the resulting equations with respect 
to the second derivatives of we obtain the following system of PDEs for the functions /*: 
fh+fa<' = 0' fk = 0' ^CC2 = 0. /cf2+/cCi = 0- g^^^^^l solution has the form = D2^+f^ 
and = —Di^ + for some smooth functions ^' = "^{t, x,y, () and /* = f^{t,x,y,(x,Cy)- 
The first summands in the expressions for /* can be neglected due to the gauge equivalence in 
the set of arbitrary elements {f^, /^). As a result, we construct the necessary representation for 
the arbitrary element H. 

Corollary 4. The class (8) is normalized. The equivalence group Gq of this class represented 
in terms of the arbitrary element H consists of the elements o/G^ with tu = and Xt = 0. The 
arbitrary elements /* are transformed in the following way: 

n'X A2 ^^^^ 

/ =>^—^-'[-x + ^x^- +y)-^)iCxc-Cy^) 

+ e{eX^f^u + na,f-^ + '-^^, 

where x = x{t) o,'^^ P = p{t^ ^) u) '^'^c arbitrary functions of their arguments. 

Proof. The class (8) is contained in the normalized subclass of the class (14) singled out by the 
constraint H(^ = 0. Therefore, any admissible transformation of the class (8) is generated by an 
element of with tu = 0, and the corresponding transformations of the space variables are 
affine with respect to these variables, Zj^ = 0. Then Dj/-' = Di{X~'^ Zf f^), i.e., the differential 
function H is a total divergence with respect to the new space variables if and only if it is a 
total divergence with respect to the old space variables. Applying the Euler operator E to the 
last equality in (15) under the conditions H(^ = 0, H^ = 0, EH = 0, EH = and tu = 0, we 
derive the additional constraint At = 0. The remaining part of the proof of normalization of the 
class (8) and its equivalence group is analogous to the end of the proof of Theorem 3. 

In order to construct the transformations of the arbitrary elements we represent the right 
hand side of the last equality in (15) as a total divergence: H = Dj/i*, where 



, nX^ 

/.^^^(/^ + A..) + ^.-(^ + f(x^ + .^))^, 

As H = Djf-' = Dj/i* = DjZ^/i*, the pair of the differential functions f^ — Zjh'^ is a null 
divergence, Dj(/-' — Z^h'') = 0. In view of Theorem 4.24 from [31] there exists a differential 
function Q depending on t, x, y and derivatives of C such that f^ — Zjh^ = — D2Q and f'^ — Zfh^ = 
T>2Q- As DjQ and, therefore, DjQ should be functions of t, x, y, C,x and C,y, the function 
Q is represented in the form Q = x(OC + p{t-,x-,y) for some smooth functions x = x{t) and 
p = p{t,x,y). □ 

Remark. The equivalence transformations associated with the parameter-functions x and p 
are identical with respect to both the independent and dependent variables, i.e., they transform 
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only arbitrary elements with no effect on the corresponding equation and, therefore, are trivial 
[22, p. 53] or gauge [36, Section 2.5] equivalence transformations. These transformations arise 
due to the special representation of the arbitrary element if as a total divergence and form a 
normal subgroup of the entire equivalence group considered in terms of the arbitrary elements 
and called the gauge equivalence group of the class (8). 

Remark. The continuous component of unity of the group Gg is singled out from Gg by the 
conditions Tt > and £ = 1. Therefore, a complete set of independent discrete transformations 
in Gg is exhausted by alternating signs either in the tuple (t, V') or in the tuple (y, -0, /"^)- 

Consider the subclass of the class (8), singled out by the further auxiliary equation /* = 0, 
i.e., the class of equations 

Ct + {i^,C} = ^if{t,Cx,Cy), C-=i^ii, (18) 

with the arbitrary elements = Cx, Cy)- 

Remark. Rewritten in the terms of H, the class (18) is a well-defined subclass of (14). It is 
singled out from the class (14) by the constraints EH = 0, = 0, Hi = and QjHqj = H. 
Indeed, the representation H = Dip{tXx,Cy) obviously implies that the arbitrary element H 
does not depend on x, y and C, is annulated by the Euler operator E and is a (homogenous) 
linear function in the totality of the derivatives dj . Hence all the above constraints are necessary. 
Conversely, the constraint EH = implies that the arbitrary element H is affine in the totality 
of Qj and, therefore, in view of the constraint CijHQ . = i/ it is a (homogenous) linear function in 
these derivatives of C,. As a result, we have the representation H = h^^Qj, where the coefficients 
h^^ , /i^^ = /i^^, depend solely on t, and Cy since H(^ = and Hi = 0. Then the constraint 
EH = is equivalent to the single equation 

whose general solutions is represented in the form h^^ = f^^, h^"^ = f^^ + f"^^ and h?'^ = f^^ 
for some differential functions /* = P{t,Cx,Cy)- This finally gives the necessary representation 
for H. 

Remark. In view of the previous remark, the subclass of the class (14), singled out by the 
constraints EH = 0, -ff^ = and i/j = is a proper superclass for the class (18) rewritten in the 
terms of H. This superclass of (18) is normalized since it is the intersection of the normalized 
class from Corollary 3 and the normalized class (8). Its equivalence group coincides with the 
group G^ described in Corollary 3. 

In a way analogous to the above proofs, the normalization of the superclass and formulas (15) 
and (17) imply the following assertion. 

Corollary 5. The class (18) is normalized. The equivalence group Gj of this class represented 
in terms of the arbitrary element H consists of the elements of G^ with tu = 0, At = 0, 
Ptt = 0, cr = and 6i = 0. The arbitrary elements /* are transformed according to (17), where 
additionally pij = 0. 

Remark. The above consideration of normalized classes is intended for the description of invari- 
ant parameterizations of the forms (8) and (18). The hierarchy of normalized classes constructed 
is, in some sense, minimal and optimal for this purpose. It can be easily extended with related 
normalized classes. For instance, the subclass singled out from the class (14) by the constraints 
EH = is normalized. Other hierarchies of normalized classes, which are related to the vorticity 
equation (5) and different from the hierarchy presented, can be constructed. 

Remark. In fact, all subclasses of generalized vorticity equations studied in this section are 
strongly normalized, cf. [36]. 
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3.4 Parameterization via direct group classification 

As proved in Section 3.3 (see Corollary 5), the class (18) is normalized. Hence the complete 
group classification for this class can be obtained within the algebraic method. Another way to 
justify the sufficiency of the algebraic method is to check the weak normalization of the class (18) 
in infinitesimal sense, i.e., the condition that the linear span the maximal Lie invariance algebras 
of equations from the class (18) is contained in the projection of its equivalence algebra (cf. 
Section 3.2) to the space of independent and dependent variables. A vector field Q in the space 
of the variables (t, x, y, tp) has the form Q = i^d^+r]d^, where the coefficients and r] smoothly 
depend on (t, x, y, -0)- For Q to be a Lie symmetry operator of an equation from the class (18), 
its coefficients should satisfy the following system of determining equations that do not involve 
the arbitrary elements (/^, /^): 

7]^^ = 0, rj^t = 0, ri^i = ^t, = "Ct , V^p - 2^1 + = 0> 

The integration of the above system immediately implies that Q £ Pg^ . 

The equivalence algebra g^ can be represented as a semidirect sum g^ = i 5 a, where 
i = >'(7^), Zix)) and a = {'Di,V2,dt,j\j^ )C{6), g{p^x + p^y)) are an ideal and a 

subalgebra of g^, respectively. Here 7^, 7^, p^, p^, 5 and x run through the set of smooth func- 
tions of the variable t and we use the notation JT^ = ^7(1), = J{t) and K,{5) = 7i{6) — Z{6). 
The intersection {kernel) of the maximal Lie invariance algebras of equations from class (18) is 

02 =('^'(7'),3^(7'),^(x))=Pi. 

In other words, the complete infinite dimensional part Pi of the projection of the equivalence 
algebra g^ to the space of variables {t,x,y,ip) is already a Lie invariance algebra for any equa- 
tion from the class (18). Therefore, any Lie symmetry extension is only feasible via (finite- 
dimensional) subalgebras of the five-dimensional solvable algebra 

a= (Pi, dt, V2, J, j')=F~a. 

In other words, for any values of the arbitrary elements /* = P{t,Cx,Cy) the maximal Lie 
invariance algebra g™^^ of the corresponding equation Cj from the class (18) is represented in 
the form g™*^^ = g^^* £ Q2, where g^^* is a subalgebra of a. A nonzero linear combination of the 
operators and ^7* is a Lie symmetry operator of the equation Cj if and only if this equation is 
invariant with respect to the algebra {J, J^). Therefore, for any extension within the class (18) 
we have that either n (J, J"*) = {0} or qJ"" D {J,J^), i.e., 

dim(gf n(J, J*)) G {0,2}. (19) 

Moreover, as Pg^ = go, the maximal Lie invariance algebra of the inviscid barotropic vorticity 
equation (5), the normalization of class (18) means that only subalgebras of go can be used to 
construct spatially independent parameterization schemes within the class (18). That is, for such 
parameterizations, the approach from [30] based on inverse group classification is quite natural 
and gives the same exhaustive result as direct group classification. Due to the normalization, 
the complete realization of preliminary group classification of equations from the class (18) is 
also equivalent to its direct group classification which can be carried out for this class with the 
algebraic method. 

Note that the class (18) possesses the nontrivial gauge equivalence algebra 
Q^--^'^ = {IC{6),g{p'x + p^y)), 
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cf. the second remark after Theorem 1. As we have PgS'^'^se — |Qj^ |^|jg projections of operators 
from gsauge obviously do not appear in g^^* for any value of /. At the same time, they are 
essential for finding all possible parameterizations that admit symmetry extensions. 

Therefore two equivalent ways for the further use of the algebraic method in this problem 
depending on subalgebras of what algebra will be classified. 

As a first impression, the optimal way is to construct a complete list of inequivalent subal- 
gebras of the Lie algebra a and then substitute basis operators of each obtained subalgebra to 
the infinitesimal invariance criterion in order to derive the associated system of equations for /* 
that should be integrated. The algebra o is finite dimensional and has the structure of a direct 
sum, = (Pi, dt, J, J^) ® (^^2)- The first summand is the four- dimensional Lie algebra gjg 
in accordance with Mubarakzyanov's classification of low-dimensional Lie algebras [27] whose 
nilradical is isomorphic to the Weyl (Bianchi II) algebra 93.1. The classification of inequivalent 
subalgebra up to the equivalence relation generated by the adjoint action of the corresponding 
Lie group on a is a quite simple problem. Moreover, the set of subalgebras to be used is reduced 
after taking into account the condition (19). At the same time, the derived systems for /* consist 
of second order partial differential equations and have to be integrated up to -equivalence. 

This is why another way is optimal. It is based on the fact that sff^ coincides with a 
subalgebra b of if and only if there exists a subalgebra b of a such that Pb = b and the 
arbitrary elements /* satisfy the equations 

en+o'fl = ^' (20) 

for any operator Q from b, where 9^ and 99* are coefficients of dt, d(^. and dji in Q, respectively. 
In fact, the system (20) is the invariant surface condition for the operator Q and the functions /* 
depending only on t and C,^ . This system is not compatible for any operator from a of the form 
Q = JC{5) + G{p^x + p^y), where at least one of the parameter-functions 6, or does 
not vanish. In other words, each operator from b should have a nonzero part belonging to 
{T>i,'D2,dt, j"^ , J^) and hence dimPb = dim b ^ 5. Taking into account also the condition (19), 
we obtain the following algorithm for classification of possible Lie symmetry extensions within 
the class (18): 

1. We classify G^-inequivalent subalgebras of the algebra d each of which satisfies the con- 
ditions dimPb = dim b and dim(b H {J,J^)) S {0,2}. Adjoint actions corresponding to 
operators from i can be neglected. 

2. We fix a subalgebra b from the list constructed in the first step. This algebra is necessarily 
finite dimensional, dim b ^ 5. We solve the system consisting of equations of the form (20), 
where the operator Q runs through a basis of b. For every solution of this system we have 
0}^* = Pb. 

3. Varying b, we get the required list of values of the arbitrary elements (/^,/^) and the 
corresponding Lie symmetry extensions. 

In order to realize the first step of the algorithm, we list the nonidentical adjoint actions 
related to basis elements of a: 



Ad(e"^*)^i = 


Vi - edt, 


Ad{e'^^)dt = 




Ad(e"^')Pi = 


--Vi+ej\ 


Ad(e"^i)J* = 


■ e-'j\ 


Ad(e"'^('^))Pi 


= Vi+elC{t5t + 5), 




= IC{e-'6{e-H)), 


Ad(e^^(^))Pi 


= Vi+eg{tpt + 2p), 


M{e'^^)g{p) 


= g{e-^'p{e-%x,y)) 


Ad(e^'^('^))9t = 


= dt + elC{6t), 


Ad(e^^*)^* = 






= dt + eQ{pt), 


Ad(e^^*)/C(5) 


= IC{6{t-e)), 
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M{e'^')dt = dt + eJ, 
kd{e'^^^^)V2 = V2 + e/C(2(5), 

Ad(e=^('^)) J* = J'+ sGitp^x - tp^y), 

where we omit tildes in the notation of operators and also omit arguments of parameter- functions 
if these arguments are not changed under the corresponding adjoint action, p = p^x + p'^y, 
ff = {p^x + p^y) cos e + {p^y — p^x) sin e, = {p^x + p^y) cos et + {p^y — p^x) sin et, 

Based upon these adjoint actions, we derive the following list of -inequivalent subalgebras 
of satisfying the above restrictions (we again omit tildes in the notation of operators): 

one- dimensional subalgebras: 

{Vi + hV2 + aj^), {dt + cV2 + cj'), {V2 + J'), {V2 + aJ^); 

two-dimensional subalgebras: 

{Vi + bV2 + aJ + lC{c) + g{cx),dt), {Vi + aj\V2 + aj^) , 
{dt + cj\V2 + aj^), {J^ + IC{6\t)),j' + IC{6\t))); 

three-dimensional subalgebras: 

{Vi+aj\dt,V2 + aJ^), {Vi + bV2,J^+IC{c\t\'^''-^),j' + ICic\t\'^^)), 
{dt + cV2,j' + IC{ce^''),j' + )C{{ct + c)e^'')), {V2, 

four- dimensional subalgebras: 

{Vi + bV2+)C{iy2),dt,J^ + }C{ui),J^ + }C{uit + uo)), {2b - = 0, buo = 0, 
{Vi,V2,j\j'), {dt,V2,j\j'); 

five-dimensional subalgebra: 

{Vi,dt,V2,j\j'). 

In the above subalgebras, due to adjoint actions we can put the following restrictions on the 
algebra parameters: a > 0, c, c G {0, 1}, a > if a = (resp. c = 0), c G {0,1} if c = 0; 
additionally, in the first two-dimensional subalgebra we can set (1 + 26) c = and ((1 + 6)^ + 
a^)c = 0; in the first four-dimensional subalgebra one non-zero parameter among uq, z/i, V2 
can be set to 1. In the last two-dimensional subalgebra, the parameters 5^ and 5"^ are arbitrary 
smooth functions of t. The subalgebras with parameter tuples ((5^, (^^) and (5^, 0^) are equivalent 
if and only if there exist constants Eq, Ei and E2 such that 5"^ = e'^'^~'^^6^{e~^^t + Eq) and 
52 = e'''6^{e-'H + Eo). 

Concerning the realization of the second step of the algorithm, we note that the system 
corresponding to the last two-dimensional subalgebra is compatible if and only if 5'^{t) = tS^{t). 
We re-denote 6^ by 5. As the general solution of the system is parameterized by functions of 
two arguments, we put the associated two-dimensional symmetry extension into Table 1, where 
the other extensions are one-dimensional. A similar remark is true for the three last three- 
dimensional subalgebras, which is why we list them in Table 2 containing symmetry extensions 
parameterized by functions of a single argument. 

The system associated with the first two-dimensional subalgebra is compatible if and only if 
(a, b) 7^ (0, —1). The solution of the system is split into three cases, (i) b 7^ —1, 1/2, (ii) 6 = 1/2 



Ad{e'''^)G{p)=g{p{t-E,x,y)), 
Ad{e'^^))C{5) = }C{e^'5{t)), 
Ad{e'^^)gip) = g{e-'p{t, e-'x, e'^y)), 

Ad{e'^')g{p) = g{n, 
Ad{e'^')g(p) = g{n, 
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and (iii) 6 = — 1 and a ^ 0. We will use the notation fj, = c/{2b — 1) for b ^ 1/2 and /x = 2c/3 
in case of 6 = 1/2. 

For the second and third three-dimensional subalgebras, the corresponding systems are com- 
patible if and only if c = c and c = 0, respectively. 

For the reason of compatibility, in the first four-dimensional subalgebra we have z^o = and 
b / —1. Due to the condition (26— l)i/i = 0, the solution of the corresponding system should be 
split into the two cases b ^ 1/2 and b = 1/2. For simplicity of the representation of the results 
in Table 3 we introduce the notation /x = z/2/ (26 — 1) if 6 7^ 1/2 and P2 = 2z^2/3 for 6 = 1/2. 



Table 1: Symmetry extensions parameterized by functions of two arguments 



-CXt 


Arguments of /i, I2 






{Vi + 6P2 + a J) 


\t\'^^^{C,xCOST + Cysinr), 
*^^nCs/COsr - CxSinr), 


r := aln \t\ 


t^"^(Ii COST - /2sinr), 
|i|*"2(Ii sinr + h cost) 


{dt + CD2 + cj') 


e''*(CxCOST + Cj/Sinr), 
e''^{C,yCOST - Cxsinr), 




e'^*(/i cos T — I2 sin r), 
e^*(/i sin r + /i cos r) 


{V2 + J') 








{V2 + aj) 






Pl,P2 


{J, J') 


t, R 




- QyP + 5{t)Qy^, 



Table 2: Symmetry extensions parameterized by functions of a single argument 



cxt 




Argument of /i, /2 




(Pi + 6p2 + aJ,5t), 




^ag(l+b)cl. 




(Pi + ip2 + a J, dt) 






R^Pl - nCy In i?, ii^-P2 + In 


{Vi-V2 + aJ,dt), 


a / 


R 


e"2*Pi - ^Cy, e-2*P2 + Kx 


{Vi + a J, V2 + aJ) 








{dt + cj\V2 + aJ) 






Pi, ^2 


{Vi+bV2,J,j') 








{dt + cD2,J,j') 






e^'\CyI^ + CxI^ - cCx<^) 


{'D2,J,J') 




t 


Pi, P2 



In Tables 1-3, Ii and I2 are arbitrary functions of two indicated arguments, arbitrary func- 
tions of one indicated argument or arbitrary constants, respectively. 

Moreover, ai = 3/(6 + 1) (for 6 7^ —1), 02 = 3/a (for 6 = — 1 and a ^ 0) and 03 = 3/(a — a) 
(for a ^ a). In Table 1, 5 is an arbitrary function of t. 



26 



Table 3: Symmetry extensions parameterized by constants 



^ext 






{Vi + aJ^,duV2 + aJ^), 


a ^ a 




{Vi + bV2,dt,J,j'),b^ 










R^Pi + {92 lnR + z^i$)Cy, R^P2 - (z>2 lnR+ ui<^)Cx 


{Vi,V2,J,j') 




t'^Pl, t-^P2 


{dt,V2,j,j') 




Pi, P2 



Up to gauge equivalence, the single parameterization admitting five-dimensional symmetry 
extension within the class (18) is the trivial parameterization, = /'^ = 0, in which we neglect 
the eddy vorticity flux. This shows the limits of applicability of the method proposed in [29], 
cf. Section 3.1. 

3.5 Parameterization via preliminary group classification 

The technique of preliminary group classification is based on classifications of extensions of the 
kernel Lie invariance algebra by operators obtained via projection of elements of the correspond- 
ing equivalence algebra to the space of independent and dependent variables [14]. It is illustrated 
here with the class (8) whose equivalence algebra Qi is calculated in Section 3.2. 

The kernel Lie invariance algebra Qi of the class (8) (i.e., the intersection of the maximal 
Lie invariance algebras of equations from the class) is {Z{x))- Denote by Qi the ideal of 
corresponding to g^, Qi = {Z{x))- In view of the classification of one-dimensional subalgebras of 
the equivalence algebra in Appendix A (list (23)) and since for preliminary group classification 
we are only concerned with extensions of the complement of in g^, we essentially have to 
consider the inequivalent subalgebras 

(Pi+aPs), {dt + hV2), {V2 + J{p)+n{a)), {J{p)+n{a)), 

{xi^') + n{a)), {n(a)+n{6) + G{p)). 

Here a G M, 6 G {—1, 0, 1}, /3 = (3{t), a = cr{t), 7^ = 7^(t) and p = p{t, x, y) are smooth functions 
of their arguments and S = 5{t,x,y) is a solution of the Laplace equation, 6xx + 5yy = 0. All 
parameters are arbitrary but fixed for a particular subalgebra. For each of the subalgebras, the 
corresponding arbitrary elements /* satisfy the equations 

efi + o'fi = ^' (21) 

where 9^ and 93* respectively are coefficients of d^, and dp in the basis element of the 
subalgebra. It now remains to present the parameterization schemes constructed, which can be 
found in the Table 4. 

In this table, Ii and I2 are arbitrary functions of four indicated arguments, r = \/ x"^ + u"^, 
If = arctany/x and Tl{a) = FTZ{a) = ^ar'^d^. 

In the last class of subalgebras no ansatz can be constructed due to the special form of 
functions Namely, as the variable ip is not included in the list of arguments of /*, any 
nonzero operator of the form TZ{a) + 7i{5) + G{p) gives an incompatible system of the form (21) 
and hence its projection does not belong to Lie invariance algebras of equations from the class (8). 

Note that some of the extensions presented are not maximal even for the general values of 
the invariant functions and I^. In particular, if an equation from the class (8) possesses a Lie 
symmetry operator of the form Af(7^) with a fixed function 7^, it possesses all the operators of 
this form. 
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Table 4: One-dimensional symmetry algebra extensions for the case /* = f^{t, x,y, (x, Cy 

Extension Arguments of Ii, I2 f^, 

{Vi + aV2) l^r^x, \t\-''v, txCx, tyCy t'^xh, r^yh 

{dt + aT>2) e-"*x, e-"*y, xCx, vCy xh, yh 



xl^ - yl^ + ^r'^Cy In r + {l3tt + at)x In r, 



{V2 + J{f3)+n{a)) 



t,ip-p\nr, ^2 

xCx + yCy, yCx - xCy 



yl^ + xP - -PCx In r + {^tt + crt)y In r 



2 



at , aC 



2/? ^ /3 



/2 + -rxy - — j-(x^ + 3xy^) 



As the class (8) is normalized (see Corollary 4), its complete group classification also can 
be obtained by the algebraic method. For this it is enough to classify only special subalgebras 
of the equivalence algebra g^", cf. a similar classification in Section 3.4 which also is used here. 
The restrictions for appropriate subalgebras are mentioned above under the classification of (at 
least) one-dimensional Lie symmetry extensions. Now we precisely formulate them: 

• The projection of a subalgebra s of gj" to the space of variables (t, x, y, V') is a Lie invari- 
ance algebra of an equation from the class (8) if and only if the corresponding system of 
equations of the form (21) for the arbitrary elements and is compatible. 

• Only subalgebras of should be classified whose projections to the space of variables 
(t, x,y,'0) are the maximal Lie invariance algebra of certain equations from the class (8). 

As a result, the classification is split into several cases. For each of the cases we have a common 
part of Lie symmetry extensions, which may be infinite dimensional. All additional extensions 
are finite dimensional and can be classified with reasonable efforts. We briefiy describe only 
the main cases arising under the classification. The complete classification will be presented 
elsewhere. 

Let s be an appropriate subalgebra of g^". As remarked above, any appropriate subalgebra 
does not contain nonzero operators of the form 'lZ{a) +Ti{6) + G{p) and includes gf' = {Z{x)) 
as a proper ideal. Denote by j the subspace of g~ spanned by the operators X{'-f^), 3^(7^), 
TZ{a), 1-L{6) and Q{p)-, where the parameters runs through the corresponding sets of functions, 
cf. Theorem 1. Then denote by rg the rank of the set of tuples of functional parameters (7^, 7^^) 
appearing in operators from 5 n j. It is obvious that rg € {0, 1,2}. We consider each of the 
possible values of tq separately. 

If tq = 0, any nonzero operator from the complement of g'^ in s has a nonzero projection to the 
subalgebra {T>i,D2,dt, , where the parameter-function f3 runs through the set of smooth 

functions of t. Suppose that the operators Q"^ = J{[i"') + with fixed linearly independent 
functions and tails T'^ = ^(7"^) + 3^(7"^) + 1z{a'') + 'H{6'') + ^(/o") G j, a = 1, . . . , n, where 
n ^ 2, belong to s. Up to G^-equivalence we can assume that = 7t((T^), i.e., 7^^ = 0, 7^^ = 0, 
6^ = and = 0. As tq = 0, we have that the commutator of any pair of operators Q's should 
be a linear combinations of certain Q's and Z{x)- This condition taken for and the other 
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Q's implies that 7"^^ = 7°^ = 0. Denote by the equation of the form (21), associated with the 
operator Q". For each o 7^ 1, we subtract the equation £^ multiphed by P"' from the equation 
multiphed by This results in the equation that does not involves /* and, therefore, is an 
identity. Splitting it with respect to Cx and Cy, we obtain the system 

i(3'(3f, - P-Pi + /3Vr - p'^aDy - P^ p% = 0, 

Taking into account that (5°^ + 5yy = and cross differentiating the two last equations of the 
system, we then derive that P^a°- — P'^a^ = 0, 5°- = 0, = 0, = and 

- P-Pl + - p-aj = 0. (22) 

Since the parameter-functions p"" are defined up to summands being arbitrary smooth functions 
of t, we can assume, in view of the equations = and Py = 0, that p"" = 0. The equations 
P^a^ — P^a^ = mean that the tuples of /3's and a's are proportional for each t, i.e., there 
exists a smooth function a = a{t) such that (cj^, . . . , cr") = a{P^, . . . , /3"). We combine the last 
condition with equations (22) and solve the resulting equations 

iP^Pf - P'^pDt + aiP^PI - P-Pl) = 

with respect to /3". The solutions are = ciaP^ f {P^)~'^adt + C2aP^ , where a = e"-^"'^* and cia 
and C2a are arbitrary constants. Therefore, the number n of linearly independent functions /3" 
cannot be greater than 2. Summing up the above consideration, we conclude that basis elements 
of 5 belonging to the complement of Qi can be assumed to have the following form: 

S^ + J{P'') + f\ b = l,...,m, JiP'')+n{a''), a = l,...,n, 

where {S^, b = 1, . . . ,m) is an m-dimensional subalgebra of {T>i,T)2,dt) and hence ^ m ^ 3, 
G j, ^ n ^ 2, the functions are hnearly independent and a"" = -(In |/3^/3| - p'^Pl\)tP'' if 
n = 2. The total dimension of extension in this case equals m + n and is not greater than 5. 

The condition ro = 1 implies that the subalgebra s contains no operators of the form J'{P)+T, 
where /3 / and T G j. Suppose that operators T' = + yij"^) + n{a') +'H{S') + Gip") 

from j, s = 1, . . . ,p, where p ^ 2 and (7''^, 7*^) are linearly independent pairs of functions, belong 
to 5. Up to Gg -equivalence we can assume that 7^^ = g, 5^ = and p^ = 0. As tq = 1, this also 
means that 7'^^ = 0, ? = 2, . . . ,p, and the parameter-functions 7'^^ = 0, s = 1, . . . ,p, are linearly 
independent. Analogously to the previous case, denote by £^ the equation of the form (21), 
associated with the operator T^. For each s 7^ 1, we subtract the equation £^ multiplied by 7^*^ 
from the equation £^ multiplied by 7^^. This results in the equation that does not involves 
/* and, therefore, is an identity. Making the same manipulations with the identity as those 
in the previous case, we obtain 6^ = 0, p'^ = 0, 7^^(t| = 7*^0"^, j^^a^ = 'y'^^a^ and, therefore, 
^l^cr^ = 7f^cr^. In view of the linear independence of 7*^ and 7^^, the last two conditions form a 
well-determined homogenous system of linear algebraic equations with respect to and cr^ and 
hence imply that a* = cr^ = 0. At the same time, if an equation from the class (8) possesses a 
Lie symmetry operator Af(7^) with a fixed function 7^, it possesses all the operators of this form. 
This means that there are only two G^-inequivalent possibility for s H j in this case, namely, 
s n j is either spanned by a single operator ^(7^^^) -|- TZ^a^), where 7*^^ and are fixed smooth 
nonvanishing functions of t, or equal to the entire set of operators of the form where 7^ 

runs through the set of smooth functions of t. Additional extensions are realized only by tuple 
of operators of the form + J{P^) + T'', b = 1, . . . , m, where G j, (5*, 6 = 1, . . . , m) is an 
m-dimensional subalgebra of {'Di,'D2, dt) and hence ^ m ^ 3. 
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Let ro = 2. We use notations of the previous case and assume summation for the repeated 
index i. Suppose that operators T*, s = l,...,p, where p ^ 3 and (7''^, 7*^) are hnearly 
independent pairs of functions, belong to s. In view of the condition rg = 2, up to permutation 
of the operators we can assume without loss of generality that 7^-'^7^^ — 7^27^^ ^ 0. Then 
for each s > 2 there exist smooth functions a'^* oft,i = 1,2, such that (7''"'^, 7'^'^) = a''* (7*^, 7*^). 
Subtracting the equation £^ multiplied by a^^ from the equation if, we derive the equation which 
should identically satisfied and, therefore, implies after certain manipulations that 5^ = a^^6^, 

= (T^ = a*V*, (t| = a*VJ and hence afV* = 0. We should separately consider two 
subcases depending on either vanishing or nonvanishing uV*. 

If (jV* 7^ then s n j coincides with the set of operators of the general form 

X{a'-i'^) + y{a'^'^) + 7^(aV*) + iiia'd') + g{a'p'), 

where (0^,0^) runs through the set of pairs of smooth functions of t satisfying the condition 
aJiT* = 0. In view of commutation relations between and operators from j, no operator of 

the form J{I3) + 7", where /3 7^ and T G j, belongs to s. Additional extensions are realized only 
by tuple of operators of the form + J{[3^) + , b = 1, . . . ,m, where S j, {S^, 6 = 1, . . . ,m) 
is an m-dimensional subalgebra of (Pi, dt) and hence ^ m ^ 3. 
Suppose that = = 0. The condition [T^,T^] G s implies that 

11 r2 I 12x2 21 f-l , 22 rl 11 2 , 12 2 21 1 , 22 1 

7 ^x+l ^y=l ^y, 7 Px+1 Py = l Px + 1 Py 

Therefore, using the push- forwards of transformations from Gq, we can set 6^ = 0, = 0. In 
other words, we can assume that the subalgebra s contains the operators = Af(7*-^) +3^(7*^), 
where 7^^7^^ — 7^^7^^ ^ 0. The system of equations of the form (21), associated with these 
operators, is equivalent to the system /*=/* = 0, i = 1, 2, which singles out the subclass (18) 
from the class (8). The complete group classification of this subclass has been carried out in 
Section 3.4. 

4 Conclusion 

In this paper we have addressed the question of symmetry-preserving parameterization schemes. 
It was demonstrated that the problem of finding invariant parameterization schemes can be 
treated as a group classification problem. In particular, the interpretation of parameterizations 
as particular elements of classes of differential equations renders it possible to use well-established 
methods of symmetry analysis for the design of general classes of closure schemes with prescribed 
symmetry properties. For parameterizations to admit selected subgroups of the maximal Lie 
invariance group of the related unaveraged differential equation, they should be expressed in 
terms of related differential invariants. These differential invariants can be computed either using 
infinitesimal methods or the method of moving frames, cf. Section 3.1. For parameterization 
ansatzes with prescribed functional dependence on the resolved quantities and no prescribed 
symmetry group, the direct group classification problem should be solved. In the case where the 
given class of differential equations is normalized (which can be checked by the computation of 
the set of admissible transformations), it is possible and convenient to carry out the classification 
using the algebraic method [36]. In the case where the class fails to be normalized (or in the 
case where it is impossible to compute the set of admissible transformations), an exhaustive 
investigation of parameterizations might be possible due to applying compatibility analysis of the 
corresponding determining equations or by combining the algebraic and compatibility methods. 
For more involved classes of differential equations at least symmetry extensions induced by 
subalgebras of the equivalence algebra can be found, i.e. preliminary group classification can 
be carried out. The framework of invariant parameterization involving methods of direct group 
classification is depicted in Figure 1. 
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Figure 1: Schematic overview of the construction of invariant discretization schemes based on 
methods of direct group classification. 
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Since the primary aim of this paper is a clear presentation of the variety of invariant parame- 
terization methods, we focused on rather simple first order local closure schemes for the classical 
barotropic vorticity equation, cf. the introduction of Section 3. That is, we parameterized al- 
ready the eddy vorticity flux v'^' using ^ and its derivatives. Admittedly, this is a quite simple 
ansatz for one of the simplest physically relevant models in geophysical fluid dynamics. On the 
other hand, it can be seen that already for this particular simple example the computations 
involved were rather elaborate. This is in particular true for the computation of the set of ad- 
missible transformations for the various classes of vorticity equations considered in Section 3.3. 
Needless to say that irrespectively of practical computational problems the same technique would 
be applicable to higher order closure schemes as well. In designing such schemes it is necessary 
to explicitly include differential equations for the first or higher order correlation terms. In the 
case of the vorticity equations, a second order closure schemes is obtainable upon retaining the 
equations governing the evolution of v'C and parameterize the higher order correlation terms 
arising in these equations. In practice, however, it becomes increasingly difficult to acquire real 
atmospheric data for such higher-order correlation quantities, which therefore makes it difficult 
to propose parameterization schemes based solely on physical considerations [43] . We argue that 
especially in such cases symmetries could provide a useful guiding principle to determine general 
classes of relevant parameterizations. 

Up to now, we have restricted ourselves to the problem of invariant local closure schemes. 
Nonlocal schemes constructed using symmetry arguments should be investigated in a subsequent 
work. This extension is crucial in order to make general methods available that can be used in 
the development of parameterization schemes for all types of physical processes in atmosphere- 
ocean dynamics. A further perspective for generalization of the present work is the design of 
parameterization schemes that preserve conservation laws. This is another aspect that is of 
major importance in practical applications. For parameterizations of conservative processes, it 
is crucial that the corresponding closed differential equation preserves energy conservation. This 
is by no means self-evident. In fact, energy conservation is violated by various classes of down- 
gradient ansatzes [45], which is straightforward to check also for parameterizations constructed 
in this paper. The construction of parameterization schemes that retain conservation laws will 
call for the classification of conservation laws in the way similar as the usual group classification. 
A main complication is that there is no restriction on the order of conservation laws for general 
systems of partial differential equations (so far, such restrictions are only known for (1 + 1)- 
dimensional evolution equations of even order and some similar classes of equations) . 

A Inequivalent one-dimensional subalgebras of the equivalence 
algebra of class (8) 

In this appendix, we classify one-dimensional subalgebras of the equivalence algebra gj" with 
basis elements (9). For this means, we subsequently present the commutator table of g^. In 
what follows we omit tildes in the notation of operators. 

Based on Table 5, it is straightforward to recover the following nontrivial adjoint actions: 



Ad(e"'^*)Pi =Vi-edt 









Ad(e^^(^))Pi =Vi+ en{t6t + 6), 
Ad{e'^^P^)Vi =Vi+ eg{tpt + 2p) 




Ad{e'^')g{p)=g{e'''p{e-%x,y)) 
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Table 5: Commutation relations for the algebra 







V2 


dt 














-dt 


Jim 


X{tll) 


©2 














-x{^') 


dt 


dt 










xijl) 


jCp) 







-jm 





-yiM + sii'Ptty) 


xm 

















3^(7') 


-3^(7?) 


-x{p^^) + g{fpttx) 


-^((7^7')*) 




-n{tdt + a-) 





-n{at) 





-H{-y'^ax)+g{-y'^aty) 




-n{t5t + 5) 


-Hixd^ + ySy - 26) 


-n{6t) 


-nifixSy - py~5^) 




Q{p) 


~g(tpt + 2p) 


-Q{xp^ + ypy+p) 




-G{Pxpy - Pyp^) 






-Z{txt + x) 


22(X) 













3^(7') 


7^(cr) 


n{S) 


oip) 


z{x) 




y{t^l) 


Ti{tat + 0-) 


n{t6t + 5) 


g{tpt + 2p) 


Z{tXt+x) 


V2 


-y{i^) 





H{x6^ + ySy - 25) 


g{xPx+yPy+ P) 


-2^(X) 


dt 


y{j?) 


7^(<T^) 






ZiXt) 


Jifi) 


x{h^) - Qi'y^Pttx) 





H0x5y - Py5,) 


g0xpy - i3ypx) 





X{j') 


^((7S')t) 


H(7Va;) - Q(ci^aty) 




g{fp.) 





yiy') 





n{fay) + Q{fatx) 


nfSy) 


gifpv) 





n{a) 


-Hij^ay) - gif^atx) 






























G{~P) 
















z{x) 


















Ad(e^^W)Pi=Pi+eZ(tXi + x), 
Ad{e'^^^'y)dt = dt + eJ{l3t), 
Ad(e^^(T'))5i = 5i + eAf(7,i), 
Ad{e'y<^^"y)dt = dt + ey{-fi), 
Ad{e''^^^'^)dt = dt + eniat), 
Adie''^^^^)dt = dt + £niSt), 
Adie'^^''^)dt = dt + eg{pt), 
Ad{e'^(^'>)dt = dt + eZ{xt), 
Ad{e'^^^''^)V2 = V2-eX{j'), 

Ad(e^^('^'))P2 = 2^2 -£3^(7'), 
Ad(e^^(''))p2 =V2 + eHixS^ + ySy - 26), 
Ad{e'^^f^)V2 =V2 + eg{xp^ + ypy + p), 
Ad{e'^^^^)V2=V2-2eZ{x), 
Ad{e'^^^'^)J{P) = A,, 
Ad(e^^(^'))J(/3) = A2, 



Ad{e'^^)Z{x) = Z{e'^x{e-'t)), 
Ad{e'^^)J{p) = jm-e)), 
Ad{e^9^)X{^^) = X{^\t-e)), 

Ad(e-9')K7')=3^(7'(t-e)), 
Ad{e'^')n{a)=n{a{t-e)), 
Ad{e^^')'H{5) = ni6it - e, x, y)), 
Ad{e'9^)g{p) = g{p{t-e,x,y)), 
Adie'9^)Zix) = Zixit-e)), 
Ad{e'^^)X{-f^) = X{e^^^), 
Ad{e'^')y{j^) = y{e'j^), 
Ad{e''^^)n{S) = 'H{e^'S{t,e-'x,e-'y)), 
Ad(e^^2)g(p) = gie-'pit,e-'x,e-'y)), 

Ad{e'^^)Z{x) = Zie'^'x), 
Ad{e'^^^^)X{j') = As, 
Ad{e'^^l^^)y{-f^) = A4, 
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Ad(e^^(^)) J(/3) = J{P) + sGiPxpy - Pyp,), 
Ad(e^^(^'));f (7^) = A'(y) + £Z((7S^)t), 

Ad(e^«(^))A'(7i) = ^'(7^) + enij'S,), 

Ad(e^^(''))-Y(7i) = X{j') + egin^p,), 

Ad(e^^('^))3^(7') = ^9, 
Ad(e^«(^))3;(72) = 37(^2) ^ e^(^2^^)^ 

Ad(e^^(^))3^(72)=3^(7^) + eg(7V,), 



Ad(e^'^(^))1W(5) = A, 
Ad(e^^(^))g(p) = ^6, 
Ad(e^'^(^'))3^(7') = y{l'') - emiW)t), 
Ad(e"'^(T'))7^((7) = ^8, 
Ad(e^'^(^'))H(<5) = H(<5(i,x - £7\y)), 

Ad(e^'^(^'))a(p) = g(p(i,x-£7\y)), 
Ad(e^^(^'))7^(a) = ^lo, 
Ad(e^y(^'))^(5) = X, y - e^^)), 

Ad(e^y(^'))g(p) = g(p(t,x,y-£7')), 



where 



^2 
^3 

^5 
^6 

^8 

^10 



J{P) + e{X{P^^) - g(/3«7'^)) + i£2z((/3(7^)2)t), 
-^(7^ cos /?£) + 3^(7^ sin Pe) - eg{j^/3u{-x sin + y cos /Se)) , 
--^(72 sin/3£) +:y(72 cos/3£) + eg (7^ /3tt(x cos /3£ + y sin ySe)), 
H((5(t,a;cos Pe + y sin/5£, — xsin/3£ + y cos Pe)), 
g{p{t, xcos /3e + y sin^£, sin/3£ + y cos /?£)), 
A'(7i) + e(7{(7W)-g(7Vy)), 
n{a) - £ (^(7Vx) - Gi^'aty)) + ^£2^ ((V )V) , 
y{j^)+e{n (7Vy) + g(7Vtx)) , 
= 7^(c7) - £ (H(7Vy) + a(7Vtx)) + ((7') V) . 



Using the above adjoint actions, we construct the following optimal list of inequivalent one- 
dimensional subalgebras of Qi'. 



{Vi+aV2), {dt + bV2), {V2 + J{p)+n{a)), ( + 7^(a) + ^(x)), 

(^-(7^) + n{a)), {n{a) + n{5) + g{p) + z{x)), 



(23) 



where a G M, 6 G {— 1) 0, 1}. In fact, each element of the above list represents a parameterized 
class of subalgebras rather than a single subalgebra. Particular subalgebras correspond to ar- 
bitrary but fixed values of parameters. Subalgebras within each of the four last classes can be 
equivalent. Thus, in the third class we can use adjoint action Ad(e'^-''i) to rescale a as well as 
the argument t of /? and a. Using Ad(e^^') allows us to shift t in the functions /3 and a. In the 
fourth class, equivalence is understood up to actions of Ad(e^-^i), Ad(e^-'^2) and Ad(e'^^*), which 
permit rescaling of cr, x and their argument t, scaling of x as well as shifts of t in /3, cr and x- 
Similar equivalence is also included in the fifth class. The last class comprises equivalence with 
respect to actions of Ad(e^'^('')), Ad(e'^'^*^''' ^) and Ad(e^'^(''' ^). In the three last classes we can 
also rescale the entire basis elements. 
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